To BE SUBMITTED TO ASTROPHYSICAL JOURNAL 
Preprint typeset using I^T^^X style cmulateapj v. 11/10/09 



DUST AND THE TYPE II-PLATEAU SUPERNOVA 2004DJ 

W. P. S. Meikle/ R. Kotak,^ D. Farrah,^ S. Mattila,"' ''' S. D. Van Dyk,'* A. C. Andersen," R. FESEN^ 
A. V. FiLippENKO,'' R. J. Foley, C. Fransson\ C. L. Gerardy,^^ P. A. Hoflich^^, P. Lundqvist'"', M. Pozzo 

J. Sollerman'^''^, and J. C. Wheeler^'' 

To be submitted to Astrophysical Journal 

Abstract 

We present mid- infrared (MIR) spectroscopy of a Type Il-plateau supernova, SN 2004dj, ob- 
tained with the Spitzer Space Telescope, spanning 106-1393 d after explosion. MIR photometry 
plus optical/near- IR observations are also reported. An early-time MIR excess is attributed to emis- 
sion from non-silicate dust formed within a cool dense shell (CDS). Most of the CDS dust condensed 
between 50 d and 165 d, reaching a mass of 0.3 x 10~^ M©. Throughout the observations much of the 
longer wavelength (> 10 /im) part of the continuum is explained as an IR echo from interstellar dust. 
The MIR excess strengthened at later times. We show that this was due to thermal emission from 
warm, non-silicate dust formed in the ejecta. Using optical/near- IR line- profiles and the MIR con- 
tinua, we show that the dust was distributed as a disk whose radius appeared to be shrinking slowly. 
The disk radius may correspond to a grain destruction zone caused by a reverse shock which also 
heated the dust. The dust-disk lay nearly face-on, had high opacities in the optical/near- IR regions, 
but remained optically thin in the MIR over much of the period studied. Assuming a uniform dust 
density, the ejecta dust mass by 996 d was (0.5 ± 0.1) x 10""* Mq, and exceeded 10~^ Mq by 1393 d. 
For a dust density rising toward the center the limit is higher. Nevertheless, this study suggests that 
the amount of freshly-synthesized dust in the SN 2004dj ejecta is consistent with that found from 
previous studies, and adds further weight to the claim that such events could not have been major 
contributors to the cosmic dust budget. 

Subject headings: circumstellar matter — dust, extinction — supernovae: general — supernovae: 
individual (SN 2004dj) 



1. INTRODUCTION 

Massive stars explode via co re collapse and eje ction 
of their surrounding layers (e.g. lArnett et al.lll989l and 
references therein). The extent to which core-collapse 
supernovae (CCSNe) are, or have been, a major source 
of dust in the Universe is of great interest. Of particular 
concern is the evidence of enormous amounts of dust 
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(> 10^ M0) in galaxies at high redshifts (z > 5). This 
comes from a variety of observations such as sub- mm and 
near-infrared (NIR) s tudies of the most di stant quasars 
(jBertoldi et all 120031: [Maiolino et aIll2004D . obscuration 
by du st of quasars in damped Ly-a systems (jPei et al.l 
|1991|) . a nd measurements o f metal abundances in these 
systems (jPettini et al.lll997l ) . Until recently, the scenario 
of dust from AGB stars tended to be rejected since it 
was thought that their progenitors woul d not yet have 
evolved off the main sequence. Howe ver, iValiante et al.l 
(|2009| ) and lDwek fc Cherchnefj (|2011l ) have argued that, 
under certain circumstances, asymptotic giant branch 
(AGB) stars may make some contribution to the dust 
budget at high redshifts. Both studies nevertheless 
cannot rule out a supernova contribution. In this paper, 
we examine the supernova option through observations 
of a nearby core-collapse event. 

CCSNe arising from short-lived Population HI stars 
might seem to be a viable alternative. It is estimated 
that each supernova (SN) must produce 0.1-1 Mq 
of dust to account for the high-re dshift observations 
(|Dwek et al.ll2?)(T7l:IMeikle et al.ll200l . Such masses have 
been predicted in mode l s of dust for mation in CCSNe 
(|Todini fc Ferraral 120011: iNozawa et al . 2003), alt hough 
more recent calculations bv ICherchneff fc Dwekl (|2010( ) 
revise such estimates downward by a factor of ~5. 
Perhaps even more problematic is the fact that actual 
dust-mass measurements in CCSNe and SN remnants 
yield values not exceeding, respectively, 10~^ Mq and 
10^^ Mq, although only a handful of such measurements 
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exist. (For a sum mary of this topic see, for example, 
iKotak et al] 120091 ^ 1.) The Spi t zer S pace Telescope 
(hereafter, Spitzer, iWerner et alj I2004D provided an 
excellent opportunity for us to test the ubiquity of dust 
condensation in a larger number of CCSNe. 

Newly-condensed dust in CCSNe can be detected by 
its attenuating effects on optical/NIR light and/or via 
thermal emission from the grains in the ejecta. Prior to 
Spitzer, the only evidence of dust condensation in typical 
CCSNe was in the Type Il-plateau (IIP) SN 1988H 
(fTuratto et al.Ml99l and SN 1999em (|Eh^iliamdi et alJ 
I2003D . However the light curve data used to type 
SN 1988H was sparse. In the case of SN 1999em, 
lElmhamdi et al.l used optical line suppression to infer 
a dust mass lower limit of about 10"'* M©. With the 
launch of Spitzer, we were at last provided with a facility 
for high-sensitivity spectroscopy and imaging of nearby 
CCSNe over the mid-infrared (MIR) range, covering 
the likely peak of the dust thermal emission spectrum. 
This can provide a superior measure of the total flux, 
temperature, and possibly dust emissivity than can 
be achieved at shorter wavelengths. Moreover, the 
longer-wavelength coverage of Spitzer allow us to detect 
cooler grains and see more deeply into dust clumps than 
was previously possible for typical nearby CCSNe. In 
this paper, we present our late-time Spitzer observations 
of the Type IIP SN 2004dj. We use these observations 
to study the dust production in this supernova. 

The paper is arranged as follows. In §1.1 we sum- 
marize and discuss previous observations of SN 2004dj. 
In §2 we present MIR (Spitzer) photometric and spec- 
troscopic observations of SN 2004dj, extending to more 
than 3 years after the explosion. This MIR coverage is 
one of the most extensive ever achieved for a SN IIP. 
We also present late-time optical and NIR photometry 
and spectra of SN 2004dj. In §3 we analyze these 
data. Corrections are derived in §3.1 for the effects of 
the line-of-sight cluster S96, and in §3.2 the mass of 
ejected ^^Ni is determined. We compare the data with 
blackbodies in §3.3, in order to assess the likely number 
and nature of the contributing sources. In §3.4 the 
origins of the IR radiation are examined in detail. The 
work is then summarized in §4. 



1.1. SN 2004dj 

SN 2004dj was discovered in the nearby spiral galaxy 
NGC 2403 on 2004 July 31 by Nakano et al. , (2004) and 
was c lassified as a normal Type IIP SN bv iPatat et al.l 
(|200l . It was the nearest such event in over three 
decades; the host galaxy lies within the M81 group. In 
IKotak et al] (|2005f ) we adopted a distance to NGC 2403 
(and the SN) of 3.13 ± 0.15 Mpc (statistical errors 
only), this being the Cep heid-derive d , zero metallicity 
value reporte d by Freedman et all (j2001| ) using the 
lUdalski et al.l (|1999( ) period-luminosity slopes. We 
contin ue to use this dis tance in the work presented 
here. iVinko et al.l ()2006D have estimated the distance 
to SN 2004dj using a combination of the Freedman et 
al. value plus their own expanding photosphere method 



(EPM) and standard candle estimates. This yields an 
average distance of 3.47 ± 0.29 Mpc, implying that the 
SN luminosity could be ~20% larger than the values 
used herein. 

The progenitor of SN 2004dj was almost certainly 
a mer nber of the compact star cluster Sandage 96 
(S96) (ISand agelll984':IB ond et al]l200l iFilippenko et al.l 
20041: iMafz -ApcUaniz et al.l 12004 iWang et alj 120051: 
Vinko et al.l I200Q) . The cluster age is various ly es- 
timated to be 14 ± 2 Myr ([Ma fz- ApcUaniz et al.l 
[200l. ~20 Mvr (iWang et al.l l2005[ ). and 10-16 Myr 
(| Vinko et al.ll2009l ). The main-seq uence mass of the pro- 
genitor is estimated at ^15 Mq (iMaiz-Apellaniz et al.l 
I200I IKotak et all I200I. ^12 M r:, (IWang et all I 2005D 
and 12 -20 M,7. (IVinko et al.ll2009D .iMafz-Apcllaniz et al.l 
(2004) and IKotak et al.l (|2005f l favor a red supergiant 
(RSG) progenitor. 



iGuenther &: Klosi (|2004D used echelle spectroscopy 
of Na I D absorption lines in NGC 2403 along the line 
of sight to SN 2004dj to infer a heliocentric velocity 
of -1-164.8 ± 0.1 kms^*. This is somewhat larger 
than the heliocentric velocity of the nuclear region of 
NGC 2403 of -130 kms'^ given in the SIMBAD and 
NED databases, but this is not surprising given the 
likely dispersion of velocities within the host galaxy. 
Indeed, [Vinko ct al. (2006) point out that H I mapping 
of NGC 2403 (jFraternah et all 12001) suggests that the 
true radial velocity of the SN 2004dj barycenter is about 
+221 kms^*. We adopt this value here. 

There is no firm consensus about the value of the 
reddening to SN 2004dj. Stellar population fitting 
for S96 yields total (Gala ctic + host) E(B - V) va l- 
ues of 0.17 ± .02 mag (iMaiz-Apellaniz et al.l 120041 ) . 
0.35 ± 0.05 m ag (IWang et al.ll2005D . and 0.1 ± 0.05 mag 
(jVinko et al.l l2009f ). Direct color comparisons of 
SN 2004dj wi th other CCSNe y ield E{B - V) values 
of_~0 mag (jZhang et al .l I2006D. 0.06 mag (|Chugail 
[200i l. and 0.07 ± 0.1 mag (| Vinko et al.l [200m. Perhaps 
most significantly, the resolved Na I D observations of 
IGuenther fc Klosd (|200l yield a host-only E{B - V) 
value of just 0.026 ± 0.002 mag. The smaller values 
of E{B — V) obtained using direct measurements 
toward SN 2004dj suggests that the SN actually lies 
ne ar the front of S96. B ased on the e x tincti on maps 
of ISchlegel et aP (fl998l ). IChugai et al] (jMl) find a 
Galactic reddening of E{B — V) = 0.062 mag while 
iZhang et all (f2006l) report E{B - V) ^ 0.04 mag (the 
same value as obtainable from NED). We therefore 
adopt a Galacti c E(B - V) = 0.0 5 ± 0.01 mag. If we 
add the IGuenther fc Klosel (|2004f l host value to the 
Galactic value, we obtain E{B - V) = 0.076 ±0.0 1 mai 
or Av = 0.24 ± 0.03 mag for a iCardelh et all (fl98^ 
extinction law with Ry = 3.1. Given the range of 
published values, in the present work we used a total 
ext inction Ay =0.31 mag, the same as that preferred 
by IVinko et al.l ()2009l ). Adoption of even the largest 
published value of E{B — V) would increase our shortest 
MIR wavelength (3.6 /im) flux by just a few percent. 

Estimates of the explosion date of SN 2004dj vary by 
several weeks. On the basis of an early-time spectrum, 
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iPatat et all ()2004[ ) placed the explosion at about 2004 
Jul y 14. This i s consi stent with the date obtained 
by iBeswick et alj (|2005[ ) , who used the radio Lp^ak 
versus rise-time relation (jWeiler et al.ll200"2l ) to yield an 
explosion date between 2004 July 11 and July 31. On 
the ot her hand, based on the EPM method, iVinko et alJ 
(|2006[) derive an explosion date as early as 2004 June 30. 
On the assumptions that the light curve of SN 1999gi 
was typical of SNe II P and that the SN_2004dj plateau 
was of similar length, iChugai et al.l (|2005l ) obtained an 
even earlier explosion date: 2004 June 13. Likewise, 
assuming si milar evolut i on be tween SN 2004dj and 
SN 1999em, iZhang et al] (|2006D find an explosion date 
of 2004 June 11. Neverthe l ess, p artly as a compromise 
with the later IPatat et al.l (|2004D value, iChugaii (|20060 
subsequent ly adopted 2004 Ju ne 28 as the date of the 
explosion. IChugai et al.l ()2007| ) also used this explosion 
date. Given the weight of evide nce for a later explo sion 
date, we reject tha t preferred by IZhang et al.l ()2006D . In 
iKotak et al.l ()2005f ) we adopted an explosion date of 2004 
July 10, or MJD = 53196.0. We use the same explosion 
date here but recognize that there is an uncertainty of 
about ±7 d. AU epochs will be with respect to MJD = 
53196.0 (i.e., < = d). 

O ptical light curv es of S N 2004d.i are p r esente d 
by iKorcakoya et al l (|2005|): 'Chugai et a l.l (|2005|): 



iLeonard et al.l (120061): IVinko et al. (2006); Z hang et al.l 
(|2006[ ): IVinko et al.l (|2009Y The VR light curves fell by 
10% and 90% of the total decline from the plateau to 
the start of the radioactive tail at, respectively, epochs 
70 ± 5 d and 96 ± 3 d. At the start of the radioactive 
tail the SN luminosity was about 15% of the plateau 
value. During the early nebular phase (up to ^300 d) 
the V band decl ined at about 1.1 mag (100 d)^^ 
(jVinko et al.ll2006[) . which is typical for a SN IIP. On 
the basis of the light curves, the mass of ej ected ^^Ni 
has been estimated at 0.02-0.03 M^d (Chugai et al.ll2Q05l: 
IKotak et al.ll2005t IVinko et al.ll200a IZhang et al.ll2006D . 
Notwithstanding, in the present work we argue that 
these authors have overestimated the ^^Ni mass and 
that the true mass is more like 0.01 M0. 

At the end of the plateau phase, SN 2004dj exhibited 
a remarkable and rapid change in some of its prominent 
optical lines, especially Ha. At 89 d the Ha profile 
still had a typical P Cygni morphology with a sym- 
metric peak blueshifted by only about —270 kms~^ 
after correction for the heliocentric velocity of the SN. 
Yet, by just ~10 d later the profile had developed a 
strong asymm etric profile with a peak at —1610 kms~^ 
(jChugai et al.l l2005i)- As time went by during the 
first year, this blueshift gradually dec reased and the 
asymmetry became less pronounced. IChugaii (|2006f l 
interpreted this unusual behavior as being due to the 
gradual emergence of an asymmetric, bipolar jet whose 
more massive component is travelling towards the 
observer. They propose that the lines are driven by 
the radioactive decay of spherical fragments of ^^Co 
cocooned in metals and he lium, lying wit h in the core. 
Using spectropolarimetry ILeonard et al.l ()2006[) also 
found evidence for departure from spherical symmetry. 
The polarization was observed to increase dramatically 



at the end of the plateau phase, implying the presence 
of significant asphericity in the inner regions of the ejecta. 

Early-time evidence of a significant circumstellar 
medium (CSM) around the progenitor of SN 2004dj 
has also been repo r ted. Radio emission was detected 
by iStockdale et al.l (|2004[) at 23 d, b y IBeswick et al. 
(|, 2005 D between 26 d and 145 d, and bv lChandra fc Ray 
(2004) at 33 d and 43 d. The SN was also detected 
at X-ray wavelengths (jPoolev fc LewinI 120041 ) at 30 d. 
The X-ray luminosity was about three times that of the 
Type IIP SN 1999em and nearly fifteen ti mes that of the 
Type IIP SN 1999gi at similar epochs. IBeswick et al.l 
(|2005D point out t hat both type s ofemission arise from 
a shocked CSM. IChugai et all ()2007D have used the 
presence of a high- velocity absorption component in the 
Ha line during the late photospheric phase to deduce 
the presence of a cool dense shell (CDS), with a mass 
of 3.2 X 10^** M0, produced by interaction of the ejecta 
with the pre-existing CSM. To this evidence for a CDS 
we add our observation of an early-time IR echo in 
SN 2004dj (see §3.4.2). 



In IKotak et al.l ()2005D we presented MIR photometric 
and spectroscopic observations of SN 2004dj at epochs 
97-137 d after explosion. Simultaneous modelling of the 
fundamental (1-0) and first overtone (2-0) of CO was 
carried out. The results favor a 15 M© RSG progenitor 
and indicate post-ex plosion CO f o rmati on in the range 
2000-4000 kms-^ Kotak et al] (|2005D also noted an 
underlying NIR continuum. A possible origin in CSM 
dust was suggested, but RSGs in S96 were favored as 
the more likely cause. Nevertheless, in the present work 
(§3.2.2) we find that the bulk of the early-time NIR and 
MIR continuum is most plausibly explained as an IR 
echo from CDS dust. The presence of 1.7 x 10"^ M© 
of Ni" *" in the ejecta was also deduced by IKotak et al.l 
(l2005h . 

In summary, while the early-time optical light curves 
and spectra of SN 2004dj are typical of a Type IIP event, 
its early radio, MIR, and X-ray behavior point to an 
exceptionally strong ejecta/CSM interaction. Moreover 
its earlier nebular-phase spectra imply an atypically 
asymmetric core. SN 2004dj is only "typical" in some 
respects. 



2. OBSERVATIONS 

2.1. Mid-Infrared Photometry 

Imaging at 3.6, 4.5, 5.8, and 8.0 /^m w as obtained 
with the Infrared Array Camera (IRAC) (jFazio et al.l 
I2004f ). at 16 and 22 /zm with the P eak-up Array (PUI ) 
of the Infrared Spectrograph (IRS) (jHouck et al.l l2004V 
and at 24 /xm with t he Multiband Imag ing Photometer 
for Spitzer (MIPS) (|Rieke et al.l [200l . Imaging ob- 
servations spanned epochs 89.1 to 1393.3 d, plus four 
observations, at 3.6 ^m only, covering 1953.9-2143.5 d. 
Aperture photometry was perforr ned on the image s 
using the Starlink package GAIA ([Draper et al.ll2003) . 
A circular aperture of radius ~ 3'.' 7 was used for the 
photometry. The background fiux was measured and 
subtracted by using a concentric sky annulus having 
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inner and outer radii of 1.5 and 2.2 times the aperture 
radius, respectively. These parameters were chosen as a 
compromise between maximizing the sampled fraction 
of source flux and minimizing the effects of the bright, 
complex background. The aperture radius corresponds 
to ~55 pc at the distance of SN 2004dj. The aperture 
was centered according to the SN WCS coordinates. 
Aperture corrections were derived from the IRAC and 
MIPS point-response function frames available from the 
Spitzer Science Center, and ranged from xl.l6 at 3.6 fim 
to x2.79 at 24 /^m. A 2a clipped mean sky estimator 
was used, and the statistical error was estimated from 
the variance within the sky annuli. Fluxing errors due 
to uncertainties in the aperture corrections are about 
±5%. 

The MIR photometry is presented in Table [TJ The 
Spitzer programs from which the imaging data were 
taken are listed in the caption. The pre-explosion MIR 
flux of S96 has not been measured, so the tabulated 
values are uncorrected for S96. An estimate of the S96 
contribution is given at the bottom of the Table. The 
effect of emission from S96 is discussed in §3.1. A tem- 
porally varying point source at the SN position is clearly 
visible in all bands. Figure [T] shows a sequence of images 
at 8.0 ^J.m at 257, 621, and 996 d. The SN is clearly 
brighter at 621 d (about x2 relative to 257 d). This is 
due to the epoch being close to the peak of the thermal 
emission from the dust (see §3.4.4.2). It can also be seen 
that the SN lies in a region of relatively bright, complex 
background emission. The MIR photometry is displayed 
as light curves and spectral energy distributions (SEDs) 
in Fig. [2] and Fig.[3l respectively. 

In Fig. [21 a rapid initial decline is seen at wavelengths 
3.6-8.0 fim. With the exception of 4.5 /im, all of the 
light curves (with sufficient temporal coverage) exhibit 
secondary maxima, with the peak emission occurring at 
epochs of ~450 d at 3.6 fim to ~850 d at 24 ^m. As 
will be discussed later, these second maxima constitute 
strong evidence of dust formation. The absence of a de- 
layed peak at 4.5 /im is due to the earlier appearance and 
dominance of CO fundamental emission in this band. Af- 
ter ~1000 d, the 24 ^m light curve starts to climb again. 
The evolution of the 24 fj,m flux is complex, as it is a 
combination of the detailed behaviors of emission from 
the ejecta dust and from an interstellar (IS) IR echo (see 
§3.4.3.1). In Fig. Owe see a steady reddening of the MIR 
SED with time. It is argued below (§3.4.4.2) that this 
effect also constitutes strong evidence of dust formation 
and cooling in the SN ejecta. The large peak at 4.5 fim 
at 114 d and 257 d is due to the aforementioned domi- 
nance of CO fundamental emission in this band. As in 
Table [U neither Fig. [2] nor Fig. [3] have been corrected for 
S96. 

2.2. Optical and Near- Infrared Photometry 

NIR imaging of SN 2004dj was obtained using LIRIS 
(Long-slit Intermediate Resolution Infrared Spectro- 
graph) on the 4.2 m William Herschel Telescope (WHT), 
La Palma, and at an effectively single epoch (spanning 
two days) with the OSU-MDM IR Imager/Spectrograph 
on the 2.4 m Hiltner Telescope of the MDM Observa- 
tory, Arizona. The wavebands are Z (1.033 yum), J 



(1.250 fim), H (1.635 /im), and K, (2.15 /zm). The 
LIRIS data were reduced using standard IRAF routinesQ 
The jittered on-source exposures were median-combined 
to form sky frames. In each band the sky-subtracted 
frames were then aligned and median-combined. 

Aperture photometry was performed oi i the reduced 
imag es using the Starlink package GAIA (jPraper et al.l 
l2002f ) with the same aperture and sky annuli as for 
the MIR photometry. The aperture was centered by 
centroiding on the sources. The magnitudes at J, 
H, and were obtained by comparison with four 
field stars lying within - 100" of SN 2004dj. The 
field-star magnit udes were acquired by measurement of 
2MASS images (|Skrutskie et al.lll997l i. For the single 
Z-band measurement, the magnitude was obtained by 
comparison with the four field stars with their JHKg 
SEDs extrapolated to the Z band. The resulting SN 
photometric measurements are listed in Table [2] and 
plotted in Fig. |4l Errors shown include uncertainties 
in the magnitudes of the four 2MASS comparison field 
stars. Pre-explosion JHKg fluxes of S96 were measured 
from the 2MASS survey (see Table[2]) and used to correct 
the JHKs light curves. Also shown for comparison is 
the 3.6 ^m light curve from the present work. In H 
and Kg the slopes flatten after ~300 d accompanied by 
a rise at 3.6 /im. This is suggestive of radiation from 
warm, newly- forming dust. 

Opt i cal p hotometry was taken from iVinko et al.l 
([20061 [2009 ). The optical light curves are displayed 
in Fig. [4l and have been corrected for emission from 
the S96 cluster using the BVRI magnitudes given 
bv IVinko et~al] (120061) . In addition, the V and R 
points of iZhang et al I (|2006[ ) around the end of the 
plateau were added to enhance the detail of this phase. 
Also shown for comparison (labelled "Rad") is the 
temporal evolution of the radioactive ener gy dep osition 
for SN 1987A as specified b v ILi et al.l (flOOl (0 to 
1200 days) and TimmeseFall (|l996[) fSOO to 3500 days) 
with the addition of the early-time contribution of ^^Ni 
decay assuming complete absorption. The radioactive 
isotopes include ^SNi, 56Co, ^^Co, ^" Co, ^ ^Na, and ^^Ti 
(In subsequent use of the ILi et al.l and iTimmes et al.l 
deposition specifications, our addition of the early-time 
contribution of ^^Ni decay is assumed.) In the optical 
hght curves, for about 35 d (115-150 d) after the end 
of the fall from the plateau the decline rate matches 
the radioactive deposition quite closely, indicating that 
this was powering the emission during this phase. The 
optical decline rates then fiatten during ~150-250 d, 
indicating the emergence of an additional power source 
(see below). After about 250 d, the optical light curves 
exhibit a steepening (possibly also present in the J band) 
which increases up to the final observations. 



2.3. Mid-Infrared Spectroscopy 

^ IRAF is written and supported by the IRAF programming 
group at the National Optical Astronomy Observatories (NOAO) 
in Tucson, Arizona, which are operated by the Association of Uni- 
versities for Research in Astronomy (AURA), Inc., under coopera- 
tive agreement with the National Science Foundation (NSF). 
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Low-resolution {R w 60-127) MIR spectroscopy 
between 5.2 and 14.5 /im was acquired at nine epochs 
between 106.3 d and 1393.3 d with the IRS in low- 
resolution mode. Long-Low (LL; 14-38 /im) observations 
were also attempted. Unfortunately, the LL observations 
were unusable. The LL slit lies at 90° to the SL slit. 
This meant that, given the scheduling constraints, the 
LL slit always lay across the host galaxy, resulting in 
heavy contamination. The MIR spectroscopic obser- 
vations were drawn from the MISC programs plus one 
epoch at 1207 d from the SEEDS program. The log of 
spectroscopic observations is given in Table [3l 

The data were processed through the Spitzer Science 
Center's pipeline software, which performs standard 
tasks such as ramp fitting and dark-current subtraction, 
and produces Basic Calibrated Data (BCD) frames. 
Starting with these data, we produced reduced spectra 
using both the SPICE and SMART v6.4 software 
packages. We first cleaned individual frames of rogue 
and otherwise "bad" pixels using the IRSCLEAN 
task. The first and last five pixels, corresponding to 
regions of reduced sensitivity on the detector, were then 
removed. The individual frames at each nod position 
were median-combined with equal weighting on each 
resolution element. Sky background was removed from 
each combined frame by subtracting the combined frame 
for the same order taken with the other nod position. 
We also experimented with background removal by 
subtracting the adjacent order. In general, nod-nod 
subtraction was preferred, as the background sampled 
in this way is expected to most closely represent the 
background underlying the SN. Any residual background 
was removed by fitting low-order polynomials to regions 
immediately adjacent to the SN position. 

One-dimensional spectra were then extracted using 
the optimal extraction tool within the SPICE software 
package, with default parameters. We found that in 
all cases the source was point-like, with a full width at 
half-maximum intensity (FWHM) that was never wider 
than the point-spread function (PSF). This procedure 
results in separate spectra for each nod and for each 
order. The spectra for each nod were inspected; features 
present in only one nod were treated as artifacts and 
removed. The two nod positions were subsequently 
combined. The nod-combined spectra were then merged 
to give the final spectrum for each epoch. Overall, we 
obtained excellent continuum matches between different 
orders. 

Despite our careful reduction procedure, the fluxes 
of the IRS spectra and the IRAC photometry were not 
completely consistent. This was due to (a) differences in 
the fixed sizes of the spectrograph aperture slits and the 
circular apertures used for the image photometry, and 
(b) the fact that the spectra were generally taken some 
days before or after the imaging data, during which 
time the SN fiux changed. We therefore recalibrated the 
IRS spectra against contemporaneous photometry in the 
8 /im band obtained by interpolation of the light curve. 
This band was chosen since it was completely spanned 
by the short-low (SL) spectrum. For each epoch, the 
IRAC 8 //m transmission function was multiplied by the 



MIR spectra and by a model spectrum of VegcQ. The 
resulting MIR spectra for the SN and for Vega were 
integrated over wavelength. The total SN spectral flux 
in the 8 /xm band was then obtained from the ratio of 
the two measurements using a zero (Vega) magnitude of 
64.1 Jy (IRAC Data Handbook, Table 5.1). These were 
then compared with the 8 /im photometry to derive 
scaling factors by which the spectra were multiplied. 
The spectra are plotted in Fig. [5j together with contem- 
porary photometric data. 

The MIR spectra comprise both continua and emis- 
sion features. Up to at least 281 d, strong emission 
from the CO fundamental was present in the IRS spectra 
and IRAC photometry. This had disappeared by 500 d. 
Strong hues of H I, [Ni I], [Ni II], [Co II], and [Ne II] 
were also present during the first year, but by 500 d only 
[Ni II] 6.64 fim, [Ni I] 7.51 /im, and [Ne II] 12.81 fim 
were still relatively strong (see Fig. [S|) . Apart from the 
CO region during 106-281 d, the 5-14 /im region was 
dominated by continuum emission. Moreover, simple ex- 
trapolation below and above the spectral coverage to the 
limits of the photometric coverage suggests that the con- 
tinuum dominated over at least 3.6-24 /im. 

2.4. Optical and Near- Infrared Spectroscopy 

We acquired optical sp ectra using ISIS o n the WHT, 
La Palma, and DEIMOS (|Faber et al.ll2003[ ) on the 10 m 
Keck 2 telescope, Hawaii. T he 895 d Keck spe ctrum 
has already been presented by iVinko et al.l ()2009D . We 
also made use of earlier post-p lateau optical spectra 
obtained by IVinko et al.l (|2006l ) a,t 89 d and 128 d 
iLeonard et all (|2006[ ) at 95 d, and iChugai et aP (|2005l ) 
at 100 d. NIR spectroscopy of SN 2004dj was obtained 
using LIRIS on the WHT and with the OSU-MDM IR 
Imager/Spectrograph on the 2.4 m Hiltner Telescope of 
the MDM Observatory, Arizona. The data were reduced 
using standard procedures in Figaro (Shortridgc 200^ 
and IRAF. The observing log for the optical and NIR 
spectra is given in Table 21 The post-100 d spectra are 
plotted in Figs. |6] and [T) 

The earlier post-plateau spectra still exhibited pro- 
nounced P Cygni features in Ha, He I 5876 A + 
Na I D, He I 10830 A, He I 20581 A, and O I 7771 A 
+ K I 7665/99 A. By 461/467 d the absorption compo- 
nents had largely vanished, with a broad-line emission 
spectrum now being observed. A few lines persisted 
to as late as the final optical spectroscopy epoch at 
925 d. We examined in detail the evolution of the more 
isolated of these lines, specifically Ha, Pa/3, [O I] 6300 A, 
[Fe II] 7155 A, and [Fe II] 12567 A. Table [5] lists the line 
luminosities (dereddened) versus epoch, toget her with 
the ra dioact ive deposition power specified by iLi et al.l 
(fT99l and iTimmes et all (Il996h for SN 1987A, but 
scaled down to an initial 0.0095 M© of ^^Ni (see §3.2). 
The evolution of the luminosities is plotted in Fig. [51 
This indicates that from just after the plateau phase to 
^460 d, the Ha and Pa/3 luminosities declined at a rate 
roughly comparable to that of the radioactive deposition. 

^ The R. K urucz Stellar Atmospheres Atlas, 1993, 
|ftp: / / ftp.stsci.edu/cdbs / grid /k93models / standards 
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In contrast, from their earliest observation at ~300 d, 
the [O I] and [Fe II] hnes decline significantly more 
slowly than the radioactive rate. Moreover, by 895 d the 
summed luminosity of just the Ha, [O I] 6300 A, and 
[Fe II] 7155 A lines exceeds that of the radioactive input 
by ~40%, rising to over 60% by 925 d. Thus, as with 
the optical light curves, we deduce the appearance of an 
additional source of energy, possibly earlier than 300 d. 

Table [5] lists profile parameters expressed as velocities 
for the more isolated lines over a range of epochs, 
shifted to the center-of-mass rest frame of SN 2004dj. 
Also listed (Col. 7) for 461-925 d are the maximum 
blue-wing velocities derived from profile model matches 
(see §3.4.4.1). Preliminary inspection indicated that 
for lines within a given element (e.g.. Ha and Pa/3), 
the velocities exhibited similar values and evolution. 
Therefore, in order to improve the temporal coverage 
and sampling. Ha and Pa/3 were grouped together, as 
were [Fe II] 7155 A and [Fe II] 12570 A. The evolution 
of the line velocities (shifted to the center-of-mass rest 
frame of the SN) is plotted in Fig. [HI and in more 
detail in Fig. [101 These plots reveal a complex velocity 
evolution. 

Up to 138 d, of the three elements considered, only 
hydrogen lines could be reliably identified. At 89 d (cor- 
responding to about half-way down the plateau-edge), 
the half width at half-maximum (HWHM) velocity was 
1820 ± 60 kms""'^, although the red and blue wings 
extended to much higher values. In addition, the peak 
emission exhibited a blueshift of —450 ± 50 kms~^. 
Then, as already described in §1, a strong asymmetry 
rapidly developed, this being attribute d to the emer - 
gence of an asymmetric, bipolar core (|Chugail [200l . 
Rather than being entirely due to the bulk motion of 
the ejecta however, some of the width of the Ha and 
Pa/3 lines may also have been produced by scattering 
from thermal electrons as in t he cases of SN 1998S 
(|Chugai 112001 and SN 2006gy (|Smith et all [2010) . but 
this effect is unlikely to have a significant influence at 
later epochs. By the time of the next observation at 
283 d, the asymmetry had diminished, with the blueshift 
of the peak now only —730 ± 50 kms~^. By this time 
the lines of [O I] and [Fe II] had emerged, also with 
asymmetric blueshifted profiles. As the SN continued 
to evolve, the line widths narrowed and by 461 d we 
see the first signs of a sharp suppression of the red 
wing. By the time of the next season's observations 
(895 d, 925 d), this suppression is very pronounced 
in all three species. This phenomenon suggests dust 
formation resulting in the obscuration of the far side of 
the ejecta and will be examined in more detail in §3.4.4.1. 

By 895 d and 925 d a second, weaker peak redshifted 
by ^170 kms~^had also appeared in the Ha and 
[O I] 6300 A profiles (Fig.[TOl). (The redshifted peak can 
also be seen in the weaker [O I] 6364 A component lying 
at -+3000 kms^Mn the [O I] 6300 A rest-frame plots.) 
The previous observations of these lines were at 467 d 
when the line luminosities were about a factor of 20 
greater (Table [5]) as well as being much wider (Table |6]). 
Therefore an underlying, weaker, redshifted component 



could have been present at 467 d or earlier, but was 
swamped by the main component of the line. The 895 d 
spectrum has a blueshift in the main peak in Ha of 
-140±10kms"i with the minor peak showing a redshift 
of about -1-160 ± 10 kms'^ In [O I] 6300 A the corre- 
sponding velocities are, respectively, — 210 ± 10 kms~^ 
and -1-170 ± 10 kms^^. Thus, the main (blueshifted) 
peak and the weaker (redshifted) peak lie roughly 
symmetrically about the local zero velocity. This 
suggests that a minor fraction of the line fiux originates 
in an emission zone centered on the SN and having 
the geometry of an expanding ring, jet, or cone. The 
[Fe II] 7155 A line also shows a secondary peak but 
at a much larger redshift of -1-480 ± 10 kms^^. There 
may actually be a peak also at around -fl70 kms~^, 
similar to those seen in Ha and [O I] 6300 A, but which 
is swamped by the stronger peak at the larger redshift. 
The profiles are analyzed in §3.4.4.1. 

We do not consider the MIR line-profile kinematics due 
to the much lower resolution (-Rfwhm ~ 100). The for- 
mation and evolution of the MIR lines will be analyzed 
in a future paper. 

3. ANALYSIS 

The evolution of the MIR spectral continuum indicates 
IR emission from dust playing a major role in the post- 
plateau flux distribution of SN 2004dj. We now make use 
of the observations described above to explore the origin, 
location, distribution, energy source(s), and nature of 
these grains. This will be done by comparison of a variety 
of simple models with the observations. 

3.1. Correction of the Supernova Flux for S96. 

The position of SN 200 4di coincides w ith that of the 
compact star cluster S96 (|Sandagdll984f) . and it seems 
likely that the progenitor was a member (|Wang et al.l 
l2005f ). Regardless of whether or not this is the case, 
it is still important to correct for the contribution of 
the cluster to the photometry and spectra, especially 
at the later epochs, before we embark on modelling the 
observed SED. To do this, in the optical region we made 
use of the pre-explosion optical photometry compiled in 
table 4 of I Vinko et all (2009). NIR photometry of S96 
was obtained from 2MASS (see Tabic [2|). Unfortunately 
there are no pre-explosion MIR images of S96 and so its 
contribution to the flux had to be assessed indirectly. 

We foimd that longward of ^0.4 /im the optical/NIR 
photometric points could be fairly represented by a 
comb ination of two blac kbodies, reddened according to 
the lCardelli et all (|1989t ) law with E{B - U) = 0.1 mag 
and Rv = 3.1. This is illustrated in Fig. [TTJ In this 
representation, the fluxes longward of ^1 /im are domi- 
nated by a component with a temperature of 3500 K. At 
shorter wavelengths, the hotter component (50,000 K) 
becomes increasingly important. This hot blackbody 
is not intended as an explanation for the shortwave 
radiation, but rather it simply serves as a means of 
representing and extrapolating the optical SED. The 
contribution of S96 to the MIR photometric points was 
then obtained by extrapolati on of the cooler bla ckbody. 
(We did not make use of the iVinko et al.l (|2009[ ) models 
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to correct for S96 as it was unclear how they should be 
extrapolated into the MIR region.) Some support for 
the effectiveness of our estimation method, at least for 
the shorter wavelengths, comes from four serendipitous 
3.6 /im images of the SN 2004dj field spanning 1954- 
2143 d (Table[T]), obtained in S'piteer program 61002 (PI 
W. Freedman). These show that, by this period, the 
light curve at this wavelength had levelled out at a mean 
value of 0.32 ± 0.02 mJy which in good agreement with 
our estimate for S96 of 0.28 ± 0.05 mJy (see Fig. [2]). 

It is possible that the above procedure could under- 
estimate a contribution from cooler material but it is 
unlikely that S96 wou l d be the source of such emission. 
iMaiz-Apellaniz et al.l (|2004f) estimate a cluster age of 
13.6 Myr and point out that by this age its parent 
molecular cloud should have been dispers e d by stellar 
winds and SN explosions. IWang et al. (I2005f) find 
an age of ~20 Myr while IVinko et al.l (|2009[ ) obtain 

10 — 16 Myr. This suggests that the flux contribution 
from S96 to longer MIR wavelengths would be small. 
On the other hand a significant SN-driven IR echo from 
the general IS dust of the host galaxy is quite likely (see 
§3.4.3.1). 

We conclude that the 3500 K blackbody extrapola- 
tion provides a reasonable estimate of the MIR flux con- 
tributions from S96. The inferred S96 fluxes in the 
3.6 — 24 ^m range are shown in Table [TJ The values are 
insensitive to the extinction over the ranges of E{B — V) 
(0.06 — 0.35 mag) suggested in the literature. It can be 
seen that the contribution of S96 at 3.6 /xm is signif- 
icant as early as ^^250 d, and dominates by ~^1000 d. 
As we move to longer wavelengths, the effect of S96 de- 
clines, becoming negligible for wavelengths longward of 
~10 /im even at the latest epochs. In the optical-NIR- 
MIR continuum modelling (see §3.4.4.2), we use the S96 
blackbody representation to correct for its contribution 
to the SED. 

3.2. Mass of Ni in the Ejecta 

It is important to establish the mass of ^^Ni in 
the ejecta of SN 2004dj since this will allow us to 
test for the presence of energy sources, other than 
radioactive decay, which might be responsible for the 
SN luminosity. In Fig. [12] we show the bolometric light 
curves (BLC s) of Vinko et al. (2006) (open circles) and 
I Zhang et al.l p006) (open triangles). The phases of 
these have been shifted to our adopted explosion date 
of MJD=531 96.0. This da te is nearly one month later 
than that of 'Zhang et al.' reducing their deriv ed ^^Ni 
mass by about 25%. The explosion date of Vin ko et al.l 
is about 10 d earlier than ours but it is not clear if this 
would significantly affect the ^^Ni mass they derived. In 
addition to the phase shifts, the BLCs of these authors 
have been scaled downward to our adopted distance 
of 3.13 Mpc. This ha s t he effect of re ducing the ^^^Ni 
masses of lVinko et al.l andlZhang et al. bv 19% and 10% 
respectively. We scaled the Vinko et al. (2006) BLC by 
a further x 1 . 1 in order to allow approximately for the 
higher total extinc tion {Ay = 0.31 mag) adopted by 
IVinko et al.l (I2009D and t he pr esent work. We scaled 
the BLC of IZhang et"all (|2006[ ) by a further factor of 
0.68 to force agreement with our adjusted version of 



the IVinko et al.l BLC. The need for this was due to 
the much larger extinction. Ay — 1.02 mag, adopted 
by IZhang et al.l . compared with the Ay = 0.31 mag 
adopted in the present work. It was found that these ad- 
justments brought our 89 d, 106 d and 129 d blackbody 
total luminosities (see §3.3 and Table [71 col. 11.) into 
fair coincidence with the other two BLCs (see Fig. IT^ . 



IVinko et all (|2006[) constructed their BLC by inte- 
grating observed fluxes in the BVRI bands and then 
extrapolating linearly from the B and / fluxes assuming 
zero flux at 3400 A and 2 3,000 A. Their BLC extended 
to 307 d. IZhang et al.l (|2006l ) simply integrated the 
observed fluxes in 12 narrow bands between 4000 A and 
10,000 A. Their BLC extended to 154 d, with flve addi- 
tional points to 180 d obtained by interpolation within a 
reduced number of bands. Thus, neither of these BLCs 
included unobserved excess flux beyond about 1 ^m. 
However, for 89-129 d the unobserved MIR flux made 
up no more than 10% of the total luminosity (see §3.4.2). 
Moreover, the optical/NIR region was dominated by 
continuum emission at this time. Consequently, the 
hot-|-warm continuum luminosities obtained from the 
present work via blackbody matching (see §3.3 and 
Table [7| col. 12.) and plotted in Fig. [T2| (solid squares) 
agree well with the two adjusted BLCs. (We exclude 
the cold component because, as argued in §3.4.3.1, 
it is due to an IS IR echo which was predominantly 
powered by the peak luminosity of the SN prior to 
the earliest epoch of observation.) By 251/281 d our 
hot-|-warm conti nuum lumino sities make up only ^60% 
of the adjusted IVinko et al.l BLC. This is due to the 
relatively strong contribution of line emission to the total 
luminosity during this time. Line emission luminos- 
ity was not included in our blackbody matches (see §3.3). 

In Fig. [121 the SN BLCs are compared with the 
rad ioactive deposit ion power in SN 19 8 7A, a s specified 
by iLi et al.l ([l993l ) & iTimmes et all ([1996!). These 
radioactive decay light curves are scaled to, respectively, 
0.0095 Mq (solid line) and 0.016 M© (dashed line) of 
^®Ni. Also shown (red) is the total radiocative luminos- 
ity in the case of 0.0095 M©. We also show (d otted lines) 
the a ctual UV-augmentcd BLCs of SN 1987A (iPun et al.l 
|1995[) derived from observations at ESQ and CTIO, 
scaled to an initial ^^Ni mass of 0.0095 M©. It can be 
seen that the 0.0095 M© case provides a good match to 
the SN 2004dj BLC during 115-150 d, just after the end 
of the plateau phase. After 150 d, unlike SN 1987A, the 
SN 2004dj BLC begins to exceed the luminosity of the 
0.0095 Mq case with the discrepancy growing steadily 
with time. Indeed, even the total radioactive luminosity 
of the scaled SN 1987A is exceeded by the SN 2004dj 
light curve, implying that an additional source of energy 
has appeared. 

Our 0.016 Mq deposition plot corresponds approxi- 
mately to the 0.02 Mq case of lVinko et all (|2006[ ) in their 
Fig. 18. We agree that this case provides a fair match 
to the BLC during ^ 260 — 310 d. Nevertheless, viewed 
within the context of the whole BLC, it can be seen 
that this "match" is actually due to an inflection section 
during the growth of the BLC excess relative to the 



8 



Meikle et al. 



true radioactive deposition. Adoption of the 0.016 M© 
case would imply an unexplained BLC deficit during 
~ 95 — 250 d. Given the phase of the event and the 
unexceptional progenitor mass it is difficult to see how 
such a discrepancy would come about. Consequently, 
we argue that only during the 115 — 150 d phase was the 
BLC of SN 2004dj actually dominated by radioactive 
decay. Beyond this period, and as deduced also in §2.2 
and §2.3, an additional luminosity source e merged. We 
therefore reject the ^^Ni mass deduced by IVinko et al.l 
(HqOi^. 

We als o reject larg e r ^^N i masses reported by other 
authors. iKotak et al.l ( 20051) used the V band exponen- 
tial tail method of lHamuvl ()2003D to derive a ^^Ni mass 
of ~ 0.022 Mq. This was based on the V magnitude at 
100 d which the subsequently more complete database 
shows was not quite yet on the radioactive tail (see 
Fig. SI , thus leading to an overestimate of t he ^^Ni mass. 
This method was also one of those used by iZhang et al.l 
dlOOe) who obtained 0.025 ± 0.010 M© of ^^Ni. As 
already indicated, their larger value was due mostly to 
their much earlier explosion epoch and m uch larger ex- 
tinctio n, neither of which we view as likely. IChugai et al.l 
(l2005l) obtained 0.020 ± 0.006 Mq of ^^Ni based on 
comparison of the V magnitude at 200 d with that of 
SN 1987A The difficul t y her e is that by this epoch (as 
also in the IVinko et al.l (20061 ) case) an additional power 
source had appeared in SN 2004dj, biasin g the derived 
^^Ni m ass to higher values. In addition, IChugai et al.l 
(|2005D used an exceptionally early explosion date, push- 
ing their re sult even high er. Finally w e note that both 
IChugai et al. (2005]) and IZhang et afl E006) al so used 
the y- light curve "steepness" method of iElmhamdi et al.l 
(j2003[) . which is insensitive to distance and extinction 
uncertainties. IChugai et al] o btained 0.013 ± 0.004 Mg, 
consistent with our result. IZhang et al.l applied the 
same method to a number of wavebands, including V, 
but obtained a larger 0.020 ± 0.002 Mq. Their steepness 
parameter at just V yields about 0.019 Mq suggesting 
that their use of multiple b ands is not th e cause of the 
apparent disagreement with IChugai et all However, the 
difference between the IChugai et al.l and IZhang et al.l 
determinations is only at the level of ~ 1.5(T significance. 

We conclude that, taking into account the uncertain- 
ties in fluxing, adopted distance, extinction and ex- 
plosion epoch (see §1.1) the mass of ^^Ni ejected by 
SN 2004dj was 0.0095 ± 0.002 Mq. We adopt this value 
for the rest of the paper. 

3.3. Comparison of Observed Continua with Blackbody 

Radiation 

Here we begin to consider the location and energy 
source of the SN continuum, especially longward of 
2 /im where thermal emission from dust would appear. 
To take an initially neutral standpoint on the interpre- 
tation, we compared optical, NIR and MIR spectra and 
photometry with blackbody continua. This provides us 
with the minimum radii of the emitting surfaces. The 
epochs were selected primarily as those for which MIR 
spectra were available, although the earliest such epoch, 
106 d, was already during the nebular era. The earliest 



MIR photometry was acquired at 89 d when the SN 
light curve was only about half way down the fall from 
the plateau to the nebular level. Given the potential 
interest of this epoch we began our model comparisons 
at this epoch despite the lack of an MIR spectrum. 
In addition, to compensate for the large gap between 
859 d and 1207 d we also considered the 996 d SED 
based on MIR photometry o nly. Optica l photometry 
was ta ken from IVinko et al.l (|2006l I2009D : IZhang et al] 
(|2006f) . Details about the sources of the other data are 
given in Tables HHH Apart from 89 d and 996 d, all the 
optical and NIR data plus the MIR photometry were 
flux-scaled by interpolation of the light curves to the 
epochs of contemporaneous MIR spectra. 

The contribution of S96, represented by a 3500 K 
blackbody of radius 1.5 x lO^'' cm (see §3.1), was first 
subtracted from all the data. To model the resulting 
0.4 — 24 /zm SN continuum it was found to be necessary 
to use three blackbodies ("hot": 5300 - 10000 K, 
"warm": 320-1750 K and "cold": -200 K). These were 
reddened and then matched visually to the continua. 
The hot blackbody was first added and adjusted to 
match the optical continuum. While the hot blackbody 
provides some information about the energy budget 
of the shorter wavelength part of the spectrum, the 
main reason for its inclusion in this study was to allow 
correction for its effect on the net continuum in the NIR 
where, up to about 500 d, it is comparable in strength 
to the warm component. The warm blackbody was 
then added and adjusted to match the — 5 — 10 ^m 
continuum plus the long wavelength end (2 — 2.4 fim) 
of the NIR continuum. It was found that, starting at 
the earliest epoch, as the SN evolved the hot-|-warm 
blackbody flux longward of 10 /im increasingly fell below 
that of the observations. Therefore a cold blackbody 
was added and adjusted to provide the final match. For 
epochs where photometry but no spectra were available 
we used the temporally nearest spectral matches to 
indicate the likely position of the underlying continuum. 
The expansion velocities of the blackbody surfaces, 
Vhot, Vwarm and Vcoid, and temperatures, Tuot, T^arm 
and Tcoid, for these matches are tabulated in Table [71 
The warm blackbody radii, Rwarm are also listed. The 
model matches are displayed in Figs. [13] (89 — 500 d) 
and [H (652 - 1393 d). In Table [7| we also show 
the luminosities, Lhot-, Lwarm and LcoM, of the three 
blackbody components together with the sum of the hot 
and warm components (L total)- We exclude the cold 
component because, as will be argued in §3.4.3.1, it is 
due to an IS IR echo which was predominantly powered 
by the peak luminosity of the SN prior to the earliest 
epoch of observation. In the final column is listed 
the radioactive deposition power corresponding to the 
ejection of 0.009 5 M© of f^Ni, s cal ed from the SN 1987A 
case specified bv lLi et al.l (|1993D fc lTimmes et al.l(|1996D . 

We stress that the blackbody matches were to the 
underlying spectral continua where this could be reason- 
ably judged, and not to the photometric points which 
also contained flux from line emission. Thus, the models 
sometimes lie below the average level of the spectra. 
This is particularly so for later epochs at wavelengths 
shortward of 2 fim where the spectra are dominated by 
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broad, blended emission lines. This tends to mask most 
of the underlying thermal continuum, leading to a pos- 
sible overestimation of the hot continuum. Nevertheless, 
the blackbody luminosities tend to underestimate the 
total luminosity as they do not allow for the total line 
emission. This is particularly the case around 200 — 500 d 
when the relative contribution of nebular line emission 
to the total luminosity is at a maximum. There is less 
of a problem before this era when the hot continuum 
is relatively strong, or afterwards when the warm/cold 
continuum increasingly dominates. Also, by 652 d, the 
relative weakening of the hot continuum means that it 
has a negligible effect at wavelengths longward of 2 ^m. 
In the 2—14 /im region the true continuum level is easier 
to judge. By 251 d and later, the total 2 — 14 fim flux ex- 
ceeds that of the continuum model by no more than 25%. 

The hot continuum declined monotonically and 
dominated the total SN continuum luminosity up to 
about a year post-explosion. It presumably arose from 
hot, optically-thick ejecta gas. At 89 d most of the hot 
continuum was probably still driven by the shock-heated 
photosphere, with the remainder being due to radioac- 
tive decay (see Table [7] and Fig.[12|). After 500 d the hot 
continuum became relatively weak and of low S/N. In 
addition, the NIR spectral dataset ended on 554 d. For 
subsequent epochs the strength of the hot component 
was estimated by extrapolation. At 500 d the hot 
blackbody match was achieved with a temperature of 
10,000 K. Between 251 d and 500 d the blackbody 
velocity declined exponentially with an e-folding time of 
~130 d and so, for epochs after 500 d the hot blackbody 
temperature was fixed at 10,000 K and the velocity 
obtained by extrapolation of the earlier exponential 
behavior (see Table [T]). While this is likely to be increas- 
ingly inaccurate with time, it is unlikely to be a serious 
source of error in determining the warm continuum; for 
example, by 500 d the hot blackbody contributed barely 
1% of the flux at 3.6 /im. Consequently, and in order 
to illustrate the MIR behavior in more detail, the op- 
tical/NIR region is not shown in the later plots (Fig.[T4|). 

The warm component luminosity declined monoton- 
ically to 281 d but then, unlike the hot component, 
increased by a factor of 3.5 by 500 d. By that time the 
hot-|-warm SN continuum luminosity exceeded that of 
the radioactive energy input by a factor of 4 and this 
excess continued to increase with time. The cause of the 
growing excess was the warm component luminosity. We 
also note that Rwarm more than doubled between 281 d 
and 500 d, but remained roughly constant thereafter. 
In fact, as we show later (§3.4.4.1, §3.4.4.2), after 500 d 
slow shrinkage in the size of the warm emission region 
occurred. 

The cold component is primarily defined by the 24 /xm 
point, and can be fairly reproduced using a range of 
temperatures (150 — 300 K) and velocities. Its luminosity 
remained roughly constant throughout the observations. 
In Table [7] we show the case with the temperature fixed 
at 200 K. 

It is interesting that, even as early as 89 d, warm 
and cold blackbodies had to be included to achieve 



a fair match to the observed fluxes in the NIR-MIR 
region. This will be discussed in §3.4.2 and 3.4.3.1. 
By 106 d, the nebular phase was just beginning and 
by 129 d the hot-|-warm luminosity was driven pre- 
dominantly by radioactive decay. By 251 d and 281 d, 
the radioactive deposition exceeded the sum of the 
hot-|-warm blackbody luminosity by about 15%. This 
excess probably went into powering the line emission 
not included in the blackbody matches. As already 
indicated, the nebular line emission was particularly 
strong at this time. Indeed it was responsible for a 
^ 50% excess of the total bolometric luminosity relative 
to the radioactive deposition (Fig. [T^ . indicating the 
emergence of an additional energy source, probably 
the reverse shock (§3.4.4.2). The appearance of an 
additional power source has already been indicated by 
the flattening of the optical light curves (§2.2). Indeed, 
after ^150 d the true BLC exhibited a steadily growing 
excess relative to the radioactive deposition (see Fig. [T^ . 

The hot and warm blackbody velocities never exceeded 
1750 kms^^, indicating continuum emission consistent 
with an origin in the ejecta or ejecta/CSM interface. The 
cold component exhibits velocities between 4500 kms^^ 
and 8500 kms~^ during the earlier phase pointing to an 
origin more likely to be outside the ejecta, specifically 
an IR echo from pre-existing dust. 

We conclude that the IR continuum comprised at least 
two components. The temperatures and temporal varia- 
tion of these components point to thermal emission from 
dust whose energy source is ultimately the supernova. 
The surge in the luminosity of the warm component by 
500 d suggests the emergence of an additional source of 
radiation — i.e., that the warm component was driven 
by different energy sources at, respectively, early and 
late times. It also raises the possibility that distinct dust 
populations were responsible for the early and late-time 
warm components. 

3.4. Origin of the Infrared Radiation 

We now explore the origin of the IR continuum radi- 
ation from SN 2004dj, especially the warm component. 
To do this we have constructed a model continuum com- 
prising hot gas, warm local dust, and cold IS dust. We 
also made use of spectral line red-wing suppression in 
the study of the warm dust. The continuum model was 
adjusted to provide visual matches to the observations. 

3.4.1. The Hot Component 

As in §3.3, the hot continuum component, presumably 
due to optically-thick ejecta gas, was represented using a 
hot blackbody having a temperature of 5300-10000 K. 
This blackbody radius and temperature was adjusted 
to obtain a match to the optical-NIR continuum. The 
warm component model (see §3.4.2 and §3.4.4.2) was 
then added and adjusted to match the ^5-10 fiia con- 
tinua plus the long- wavelength end (2-2.4 fim) of the NIR 
continuum. As with the pure blackbody matches it was 
found that, at all epochs, the warm component model 
flux longward of 10 ^m increasingly fell below that of 
the observations. This excess is attributed to an IS IR 
echo (see §3.4.3.1). 



10 



Meikle et al. 



3.4.2. The Warm Component: Early Phase IR Excess due 

to a CDS 

As pointed out above, a striking result from the 
hot-blackbody matches is that an NIR-MIR excess 
(relative to the hot component flux) was present as early 
as 89 d post-explosion. (For brevity we henceforth refer 
to the NIR-MIR continuum excess as the "IR excess".) 
The MIR spectra from 106 d onward show that this was 
primarily due to continuum emission. We have no reason 
to suspect that the IR excess in the 89 d SED was not 
also due to continuum emission. Indeed, the CO peak at 
4.5 /im is noticeably suppressed at this epoch, compared 
with 106 d (Fig. [3|). The obvious interpretation of the 
IR excess is that it arose from warm dust heated by the 
supernova. The very early appearance of this emission, 
when the H-recombination front had not yet reached the 
He-metal core, argues against an origin in newly-formed 
ejecta dust. In addition, SN-ei ecta dust-formation 
models fe.g. lTodini fc Ferraral[200ll iNozawa et al.ll2003[ ) 
suggest that dust formation in a CCSN is unlikely to 
occur until after one year post-explosion. A second 
possibility of direct shock-heating of pre-existing cir- 
cumstellar dust is also ruled out. The SN UV flash 
would evaporate dust out to (0.5 — 1.0) x 10^^ cm 
(0.016 - 0.032 pc) (see §3.4.3.1). To reach this distance 
by 89 d would require a shock velocity of at least 
65,000 kms~^. Yet the velocities of th e Ha trough 
(jKorcakova et al.l 120051 : iVinko et al.]|2006[ ) as well as of 
the extren ie blue-profile edge (our measurement s of the 
spectra of iKorcakova et al.i (|2005l ): IVinko et al.l ([2006)) 
during the period 25 — 90 d suggests that the bulk of 
the ejecta never exceeded velocities mu c h mo re than 
~15,000 kms-^ Indeed iChugai et all (|2007[ ) adopt 
a velocity of 13,000 kms~^as the boundary velocity 
in their treatment of SN 2004dj. A third possibility, 
CSM dust heating by X-rays, is also implausible. The 
luminosity of the warm continuum component at 89 d 
was 1.3 X 10''° erg s~^ (Table [7]) but the X-ray luminosity 
at 30 d was only about 1% of this (|Poolev fc LewinI 
[200l . In the case of dust in a CDS, collisional heating 
is ruled out as the energy available in the CDS is about 
a factor of 10^ less than that required to account for the 
observed IR excess. The heat capacity of the grains is 
also insufficient to account for required energy. 

The most likely explanation for the early-phase IR 
excess is an IR echo of the SN early-time luminosity 
from circumstellar dust. There are two possible sce- 
narios here. In the first of these, the IR echo is from 
pre-existing dust in the CSM. Such early IR echoes have 
been suspected before in other CCSNe. IWood en et al.l 
([l993) found an IR excess in SN 1987A as early as 
260 d and possibly also at just 60 d. They hypothesised 
that the origi n of the e xcess was warm, SN-heated dust 
in the CSM. iWoodenI ()1997f ) reiterates th at the cause 
was C SM dust "echoing the light curve". iFassia et al.l 
(|2000( ) reported a strong K — L' excess in the emission 
from the Type Iln SN 1998S at 130 d. They attributed 
this to pre-existing CSM dust heated either by the SN 
luminosity (a conventional I R echo) or by X-ra ys from 
the CSM-shock interaction. iPozzo et al.l (|2004l) argued 
in favor of the former scenario. 



In the second IR echo scenario, as the fast moving 
ejecta collides with the CSM, a CDS form s between the 
forward and reverse shocks. As noted in ?1. IChugai et al.l 
([2007) have used the Ha spectrum to deduce the exis- 
tence of a CDS in SN 2004dj. ^ yithin the CDS con ditions 
can allow new dust to form (jPozzo ei~an[200l . CDS 
dust has been invoke d as the origin of the early-time IR 
exces s of SN 2006ic (IS^ith et al.l [20081 iDi Carlo et al.l 



l2558HMattila et al.i[200'8[ ). iMattila eTahl (|2008f ) showed 
that the IR emission was probably an IR echo from the 
CDS dust. IR excesses in CCSNe at later times have 
also been attributed to emission from newl y- condensed 
CDS dust viz. in the Type Iln SN 199 8S dPozzo et al " 
20041) and the Type IIP SN 2007od ([Andrews et al 



2010 ) 



We explored the possibility that the early IR excess 
was due to an IR echo of the SN luminosity from 
spherical distributions of either pre-existing circum- 
stellar dust or newly- fo rmed CDS dust. Details of 
the model are given in iMeikle et al] ()2006t ). It fully 
allows for the effects of light-travel time across the 
dust distrib ution. Vers ions of this model have als o 
been used in IMeikle et a l. (2007); Mattil a et all ([2008); 
iBotticella et al.l (|2009D : iKotak et al. (20091 ). The model 
assumes a spherically symmetric cloud of grains cen- 
tered on the SN, with a concentric dust-free cavity 
at the center. The SN is treated as a point source. 
For simplicity, a single grain radius, a, is adopted. 
For ease of computation, we assumed that the grain 
material was amorphous carbon where, for wavelengths 
longer than 27ra, the grain absorptivity/emissivity can 
be well-approximated as being proportional to A~^'^^ 
(jRouleau fc Martin Ill99l() . For shorter wavelengths, an 
absorptivity/ emissivity of unity was used. The m aterial 
density is 1.85 gcm~^ (|Rouleau fc Martini 119911 ). Free 
parameters are the grain size, grain number density, 
radial density law and extent. The input luminosity is a 
parametrized description up to 550 d of the BLC shown 
in Fig. 112! viz.: 

Lho\ ~ Loexp{—t/T), where 
Lo = 57.0, T = 1000.0 d for t < 0.2 d, 
Lq ^ 1.70, T ^ 23.4 d for 0.2 <t< 23.0 d, 
Lq ^ 0.66, T ^ 141.9 d for 23.0 < t < 76.4 d, 
Lo = 50.7, T = 15.9 d for 76.4 < t < 112.0 d. 
Lq ^ 0.079, r = 174.8 d for 112.0 < t < 550.0 d 
Lo = for t > 550.0 d 

Lq is in units of lO''^ ergs^^. The brief but highly 
luminous first term represents the energy in the UV 
flash. The second term, which covers the pre-discovery 
phase, was estiniated b y using the BLC of SN 1999em 
(jElmhamdi et al.l I2003D . adjusted so that the epoch 
and luminosity of the beginning of its plateau phase 
coincided approximately with the earliest observed point 
on the SN 2004dj BLC. 

For the third, fourt h and fifth terms, the adjusted 
BLCs of IVinko et al.l (l2006l ) and IZhang et al.l (|2006l ) 
were used (see §3.2). As already pointed out, these 
BLCs did not include excess fiux beyond about 1 /xm, 
nor any flux beyond 2.3 /im; virtually all of the warm 
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and cold components were excluded. Thus, use of 
the above parametrized description avoids "double- 
counting" of a putative CSM/CDS IR echo. Indeed the 
BLC description slightly underestimates the true SN 
luminosity input as it excludes all IR emission longward 
of 2.3 fj,m, as well as all UV emission shortward of 
0.34 /im. Consequently the input luminosity was scaled 
by a factor of about 1.1 to allow for unobserved UV and 
IR fluxes. 

For the case of pre-existing CSM dust the model 
was adjusted to reproduce the IR-excess SEDs for the 
three earliest epochs (89, 196 and 129 d). The outer 
limit of the circumstellar dust was initially set at 10 
times that of the cavity radius, and a (steady 
wind) density profile was assumed. However, it was 
impossible to reproduce the quite rapid temporal decline 
of the IR excess without raising the temperature of 
the hottest grains to above their evaporation temper- 
ature. While a better match to the SED shape and 
evolution was obtained by setting the density profile to 
steeper than r~^, the best match was achieved with a 
discrete, thin shell. We therefore adopted this configura- 
tion, setting the shell thickness at xO.l the cavity radius. 

We investigated a range of cavity radii. For pre- 
existing CSM dust the minimum size of the concentric 
dust-free cavity is fixed by the extent to which the dust 
was evaporated by the initial UV flash from the super- 
nova. In this scenario, while the bolometric light curve 
(BLC) dominates the heating of the surviving dust, the 
size of the dust-free cavity is determined by the lumi- 
nosity peak of the UV flash, with Vevap oc L^f^j. (|Dwekl 
119831). For a Type IIP SN the flash luminosity is esti- 
mated to peak at about 10 "^^ ergs~^ ()Klein fc Chevalieii 
Il978t iTominaga et al.l[2009( ) although it has never been 
observed directly. A similar peak luminosity is estimated 
for the Type Ilpec SN 1987A (Ensman & Burrows 1992). 
iDwekl ()1983[ ) provides an approximate estimate of the 
flash-evaporated cavity size for a Type II supernova . 
More recently, in a detailed studv iFTschera et "all ()2002f ) 
determined that for SN 1987A the UV flash would have 
totally evaporated graphite CSM dust out to a radius of 
(0.5 - 0.9) X 10^'^ cm (0. 016 - 0.029 pc). We therefore 
adopt the iFischera et al.l estimates and apply them to 
the case of SN 2004dj. In any case, matches to the data 
are fairly insensitive to the cavity radii. Fair matches 
to the early-time IR-excess SEDs were obtainable with 
cavity radii (0.5 - 1.3) x 10^^ cm (0.016 - 0.042 pc) and 
corresponding grain radii of 0.1 — 0.04 /zm. For cavity 
radii exceeding ~0.05 pc the model continmmi slopes 
were inconsistent with the observations. 

As an example of a pre-existing CSM dust model 
match we have: 

cavity radius: 0.7 x 10^^ cm (0.023 pc), 
shell thickness: 0.7 x 10^^ cm (0.0023 pc), 
grain size: 0.07 /xm, 

grain number density: 6.0 x 10"^ cm~'^, 
total grain mass: 0.38 x 10~^ Mq. 

The optical depth through the shell in the optical band 



is 0.006, which is easily encompassed within the observed 
total extinction of Ay = 0.3 mag. For a dust/gas mass 
ratio of 0.005 (see §3.4.3.1), the dust mass corresponds 
to a total shell mass of 0.76 x 10~^ Mq. For the adopted 
shell thickness and a typical RSG wind velocity of 
20 kms^^, this mass would be produced by a mass loss 
event about 1100 years ago, lasting for 110 years, with 
a mass loss rate rate of 7 x 10~^ Mq yr~^. However, 
while this mass loss and loss rate are plausible for 
RSGs, such discrete events are not thought to occur 
in this type of star. This prompts us to seek a more 
natural explanation for the thin dust shell. An obvio us 
candidate is the CDS inferred bv lChugai et al.l (|2007| ). 

For the CDS scenario the size of the cavity i s essentially 
the radius of the CDS. From .Chugai et al.l f2007.) . the 
CDS radius is given by 

rcDS = 5.2 X 10^^ X {{t + 11.5)/64.0)^/^ cm (1) 

where t (in days) is with respect to our explosion epoch. 
Thus, at the earliest of our epochs, 89 d, the CDS 
radius was just 7.7 x 10^^ cm, or about one tenth of 
the pre-existing CSM dust cavity. This is only 3 light 
days implying that light travel time effects are small. 
Nevertheless, for convenience and ease of comparison 
with the pre-existing dust case, we applie d the IR echo 
model to the CDS case. Although the iChugai et al.l 
study terminates at just 99 d, we assume that the CDS 
radius continues to increase as described in equation (1) 
until at least 500 d, by which time its contribution to 
the IR excess is small. 

We compared the CDS case of the IR echo model with 
the observations, with the dust lying in a thin shell of 
radius reus as given above. We found that setting the 
dust mass at a constant value produced a poor match to 
the observations. With a match at 106 d, the model at 
89 d yielded a continuum which matched the observed 
IR excess at 8 /im but exceeded the IR excess by nearly 
x2 at 3.6 /im. On the other hand, at 129 days, while 
the model continuum slope was similar to that of the 
observed IR excess, it significantly underproduced it 
with the deficit being as large as ~35% at 8 /^m. We 
propose that this problem is due to the unjustified 
assumption of a fixed CDS dust mass. CDS dust could 
not form until the supernova flux at the CDS had faded 
sufficiently for proto-dust material to cool below the 
condensation temperature. This occurred at about 50 d, 
assuming amorphou s carbon dust. Therefore, following 
iMattila et"all ()2008D we allowed the CDS dust mass to 
grow as Md = Mo{l — exp{—{t — to)/td)) where t is time, 
to is the time at which dust condensation began (set 
at 50 d), td is the characteristic grain growth timescale 
and Mq is the asymptoticall y-approached flnal mass. 
(We note that IMattila et al.l also deduced an epoch 
of 50 d for the start of the CDS dust condensation 
in SN 2006jc.) No attempt was made to simulate the 
growth of individual grains which were assumed to 
appear instantaneously at their flnal size. Owing to the 
light travel time differences across the CDS, the grain 
condensation is seen to commence during the epochs 
(^0 - {rcDs/c)) to {to + (rcDs/c)) days. Yet, even as 
late as 500 d, {rcDs/c) was only about 12 light days 
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and so this effect was ignored. 

In Fig. [151 compare the CDS model hght curves 
and spectra with, respectively, the observed MIR excess 
fluxes at 3.6 /im and 8.0 /im (LH panel), and with 
the 89, 106 and 129 d MIR excess SEDs (RH panel). 
The free parameters for the CDS case were the dust 
mass scaler, the grain radius and the grain growth 
timescale. Satisfactory reproduction of the IR excess 
was achieved for all five epochs spanning 89-281 d with 
td = 50 d, a grain radius of 0.2 /im and a final dust 
mass of Mq = 0.33 x 10~^ M0, similar to the dust mass 
derived in the pre-exist ing CSM case abov e. The CDS 
mass is 3.2 x lO"'^ Mq (|Chugai et al.ll2007[ ) indicating a 
plausible final dust/gas mass ratio of 0.01. The hottest 
dust ranged from 1330 K at 89 days declining to about 
640 K at 281 days. The optical depth in the UV/optical 
range was 0.077 at 89 d faUing to 0.011 by 500 d — 
i.e., consistent with the observed total extinction of 
Ay — 0.3 mag. By 500 d the MIR excess exceeds that 
of the model by a significant factor (especially at longer 
wavelengths) (Fig. I15p implying the appearance of an 
additional energy source. The CDS contribution to the 
SN continua for 89-500 d is plotted in Fig. [HI The 
luminosity contribution of the CDS IR echo, Lcds^ is 
listed in Table 9, Col. 3. For epochs after 500 d the CDS 
component was negligible and so was ignored. 

We conclude that the early-time IR excess was proba- 
bly due primarily to an IR echo from newly-formed dust 
lying within the CDS. The rapid decline of the IR excess 
flux, especially between 89 d and 106 d, is due largely to 
the ongoing fall from the BLC plateau, tempered by the 
growth of CDS dust during this period. Owing to the 
small size of the CDS, light travel time effects are small. 
In contrast, in the pre-existing CSM dust scenario, the 
rapid decline of the input luminosity is tempered by the 
larger size of the dusty shell which produces significant 
light travel time effects, smoothing the observed IR 
excess light curve over longer tiniescales. Nevertheless, 
given the natural explanation by the CDS model of the 
required thin shell dust distribution, in the completion 
of the analysis below we use the CDS scenario. 



3.4.3. The Cold Component 
3.4.3.1. An Interstellar Echo 

The cold component is defined primarily by the 24 /^m 
data and can be fairly reproduced using a range of 
temperatures (150 — 300 K) and velocities. As described 
above (§3.3), the luminosity of the cold component 
remained roughly constant throughout the observations. 
In addition, the cold blackbody velocities were high, 
arguing against an origin in ejecta dust. An obvious 
alternative is an IS IR echo. The possibility of detecting 
the refiection of S N optical light from IS grains was 
first suggested by Ivan den Berghl ([19651 ). In the IR 
a potentially much more important phenomenon is 
the absorption and re-radiation by the grains of the 
SN BLC energy (the "IR echo"). The possibility of 
detecting an IS IR echo from a SN was first proposed 



bv IBode fc E"to^ (fT980l) l^ The occurrence of cold dust 
IS IR echoes should be relatively common for SNe 
occurring in dusty, late-type galaxies. The SEDs of such 
echoes tend to peak in the 20 — 100 /im region, allowing 
echo detection in nearby galaxies by Spitzer during its 
cold mission. Spitzer-hased evidence of this phenomenon 
in the Milky Way Galaxy have been presented for the 
Cassiopeia A SIV (iKrause et al.l [200I iKim et all 120081: 
iDwek fc Arelldtl I2008D . iMeikle et al.l (|2007D showed 
that an IS IR echo provided a natural explanation for 
the strength and decline of the 24 /im flux between 
670 - 681 d and 1264 d i n SN 2003gd in th e SA(s)c 
galaxy NGC 628 (M74). IKotak et al.l ((2009!) showed 
that the cold component of the SN 2004et SED was most 
likely due to an IS IR echo in the SAB(rs)cd host galaxy 
NGC 6946. The host galaxy of SN 2004dj, NGC 2403, 
is of type SAB(s)cd and so there is a good likelihood of 
a similar IS IR echo occurring. Therefore, we included 
an IS IR echo component in our modelling of the SED. 

Our IS IR echo model is the same as that used for 
the early-time CDS IR echo (§3.4.2). Only the dust 
distribution and grain radii are different. We ignore 
the CDS dust emission derived in §3.4.2 since this is 
already invoked in the CDS echo model and included 
in the continuum modelling up to 500 d. We recognise 
that a more extended, lower density CSM may also have 
existed. However, we found that the addition of yet 
another model component was unnecessary to provide 
plausible matches to the continua and so for simplic- 
ity the possible effects of an extended CSM were ignored. 

As described in §3.4.2, pre-existing dust surrounding 
the SN would have been evaporated by the SN flash 
out to a distance of about 0.025 pc. At sufficiently 
late epochs, enlargement of this cavity can be pro- 
duced by the forward shock. Assuming that the shock 
velocity is comparable to the highest ejecta velocity 
viz. ~15,000 kms~^ (see §3.4.2), and that there was 
no deceleration, the edge of the UV-fiash-determined 
cavity would be reached after about 600 d. After 
this time, the SN shock would evaporate the dust and 
enlarge the dust-free cavity. Therefore, for epochs 
earlier than 600 d we fixed the inner limit of the IS 
dust at 0.025 pc. For later epochs the inner limit of 
the IS dust increases from 0.028 pc (34 light days) at 
662 d to 0.058 pc (70 light days) by 1393 d. It might 
be objected that this ignores the possibility of shock 
deceleration. However, the IR echo contribution which 
the shock-evaporated dust would otherwise have made, 
for 600 > t > 1393 d, to the cold component is neghgible 
since the cavity radius was never more than 70 light 
days i.e. during this late period the SN peak luminosity 
would have long since passed by the dust near the cavity. 

The outer limit of the IS dust is much less cer- 
tain. iBendo et al.l ([2O1O0 have used Spitzer imaging 
of NGC 2403 at 70 /im and 160 /im to map the dust 
column density via its thermal emission. They also used 
the far-IR images in conjunction with H I observations 
to map the gas/ dust ratio. The objection to using such 

^ IWrigh3 1)19801 ) also considered an IR echo from a SN but only in 
the more restricted case of an explosion within a molecular cloud. 
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measurements in the present work is that they only pro- 
vide total densities through the disk (or to the mid-plane 
assuming symmetry), and not directly to SN 2004dj. 
While it was argued above that SN 2004dj actually lies 
towards the front of S96, we are still faced with the 
uncertainty of the depth of S96 in the galaxy plane. 
Our approach, therefore, is to assume a spherically 
symmetric dust distribution centered on the SN, with an 
outer limit of 100 pc. While not appropriate in principle 
for the outer limits of the dust in the galactic disk, 
for the early era being considered spherical symmetry 
provides a good approximation; at this epoch, the echo 
ellipsoid is extremely elongated and so the region of 
the spherical model outer surface intercepted by the 
ellipsoid is small. The plane of NGC 2403 appears to be 
tilted at an inclination of about 55° roughly doubling 
the face-on column density, and so the adopted outer 
limit is equivalent to a ~60 pc scale height for the IS 
dust with the host galaxy face-on. 

The free parameters were (i) the grain size, which 
influenced the dust temperatures, and (ii) the grain 
number density, which determined the luminosity. 
These parameters were adjusted in conjunction with the 
warm dust models (see §3.4.2 and §3.4.4.2) to provide 
a match to the longwave excess. It was found that 
satisfactory matches at all epochs were obtained with 
a grain radius of 0.1 /zm and a dust number density 
of 2.6 X lO"^'^ cm"-^ (i.e., an IS gas number density 
of 0.24 cm"-^ and a dust/gas mass ratio of 0.005 - 
iBendo et al.ll2010l ) This is comparable to the typical 
value of 5 X 10" cm'^ for the Milky Way jMi^Ell) 



and 7.0 x 10"" cm"^ obtained bv iKotak eFall (|2009l ) 
for NGC 6946. In the wavelength region (A > 15 /im) 
where the cold component makes a significant contribu- 
tion (> 20%) to the total flux, the IS IR echo model 
(15-150 ^m) maintained a near-constant luminosity of 
~ 2.0 X 10^^ ergs-^between 89 d and 1393 d. 



The IBendo et al.l dust column density map of 
NGC 2403 indicates 0.05 Mgpc'^ at the position of 
SN 2004dj. From the above dust number density from 
the echo model and integrated over 100 pc we obtain 
jus t 0.003 M(np c~^. The absorption opacities used 
by IBendo et al.l are consistent, to within a factor of 
~ 2, with the opacity law used in the present work. 
Thus, w e have a ^ x 8 discrepancy in column density 
between IBendo et al.l and the present work. Part of 
this discrepancy may be that the adopted dust outer 
limit is too small. The derived IS dust density is fairly 
insensitive to the extent of the outer limit. For example, 
with a 75% increase in the outer limit the model match 
to the longwave continuum is retained with a reduction 
of just ~10% in the IS dust density. Thus, at least 
some of the discrepancy could be removed by simply 
increasing the dust outer limit. An explanation for the 
remaining discrepancy is that the SN and S96 actually 
lie significantly above the mid-plane of NGC 2403. 
This is confirmed as follows. The optical depth to 
UV/optical photons yielded by the model is about 
0.026. We can use this to estimate independently 
the host extinction to SN 2004dj. The optical depth 
translates to an absorption-only extinction of 0.028 



magn itudes. Assuming an albedo of about 0.4 (|Drainel 
l200l . we obtain a total host extinction of 0.071 mag. 
This is reasonably consistent with the host-only value of 
Ay = 0.081 mag ( E = 3.1) obtained v ia the Na I D ob- 
servations of iGuenther fc Klosd ()2004l ). A s imilar result 
is ob tainable from our study of SN 2004et (jKotak et al.l 
l200a . For SN 2004et, the IS IR echo model yields 
an absorption-only extinction of 0.081 mag implying a 
total host extinction of about Ay = 0.20 mag. High 
resolution spectra of Na I D lines to SN 2004et gave 
a total (host plus Galaxy) E{B — V) = 0.41 mag 
(jZwitter et al.ll2004f ). or Ay = 1.27 mag (Rv = 3.1). 
Subtracting the estimated Galactic contribution o f 
Ay = 1.06 mag (jSchlegel et al.lll998t iMisra et al.ll2007D 
yields a host-only value of Ay = 0.21 mag, in excellent 
agreement with the IS IR echo-based value. 

One possible objection to the IS IR echo interpretation 
of the cold component of SN 2004dj is that the steady 
component of the flux around 24 /xm could be due to 
cool IS dust in S96, insufficiently corrected for by the 
3500 K blackbody extrapolation ( §3.1). W c regar d this 
as unlikely. As already indicated. iMaiz-A pcllaniz et al.l 
pool ) argue that the parent molecular cloud of S96 
should have been dispersed by stellar winds and SN 
explosions. In addition we have found that (a) a single 
set of IS IR echo parameters provided a fair match 
throughout the 89 — 1393 d covered and (b) the derived 
dust column density is comparable to that obtained 
from IS Na I D spectroscopy. 

A second possible objection to the IS IR echo inter- 
pretation is that the cold component is actually due to 
free-free emission. This is discussed and dismissed in 
§3.4.3.2. 

We conclude that the cold component of SN 2004dj 
was due to an IS IR echo. As originally suggested by 
iBode fc Eva!ni (|1980f) . the study of IS IR echoes from 
SNe can provide an independent method of measuring 
the IS extinction in galaxies. 



3.4.3.2. Free-Free and Free-Bound Radiation 



In their MIR study of SN 1987A, iWooden et all ([l993l ) 
suggest that during the 60 — 415 d period, a significant 
proportion of the MIR flux longward of ~20 /zm was 
due to free-free (ff) emission, with ff and free-bound (fb) 
emission also contributing at shorter wavelengths. For 
SN 2004dj we have looked at the possible contribution of 
fb-ff radiation. It is particularly important to consider ff 
emission as this rises in flux towards longer wavelengths, 
thus potentially reducing or even dismissing an IS IR 
echo contribution. For consistency in this analysis fb 
emission must also be included. For the wavelengths 
covered this is strongest in the optical/NIR region. 
Free-bound radiation has a strong "sawtooth" structure 
and the extent to which this structure is undetectable in 
the optical/IR continuum places a limit on its strength. 
For a particular temperature this, in turn, constrains 
the strength of the ff emission {Fff/Fft, cx ^/T/v^). 
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The fb-ff radiation was assumed to arise in a hydrogen 
envelo pe centered on the supernova. iKozma fc FranssonI 
(|1998D modelled the time-dependent behavior of the tem- 
perature and ionization of SN 1987A, including the hy- 
drogen envelope, for epochs after 200 d. Therefore, we 
made use of their study to model the ff and fb emission 
from SN 2004dj at epoch 251 d and later. We adopted a 
hydrogen envelope of a few solar masses and used their 
density profile: 

p = 9.1 X 10"i^(t/(500 d))-^{v/{2000 kms"^)"^ g ^^-3 

(2) 

By 251 d, the highest velocity of observable hydrogen in 
SN 2004dj was no more than - 3100 km s^^' (Table i]). 
By 925 d the maximum H velocity observed in SN 2004dj 
was no more than 1400 kms~^. There were no optical 
or NIR spectra covering the final three MIR epochs. 
We therefore assumed that the maximum hydrogen 
velocity at these late epochs remained at 1400 kms~^. 
The hydrogen inner limit is less certain; a value of 
1000 kms~^was adopted for all epochs. In any event, 
the fb-ff luminosity from 500 d onwards was completely 
negligible for any plausible hydrogen velocity limits (see 
below) . 

We estimated the fb-ff emission from the hydrogen 
within the velocity limits. This was don e using the 
escape probability formalism (Os terbrocld !l989), al- 
though in practice this was unnecessary as the matches 
to the data always showed that the hydrogen was very 
optically thin. The mass of the fb-ff-emitting hydrogen 
feh from ~4 Mq at 251 d to 0.75 Mq at 996, 1207 
and 1393 d. Using i nodes t extrapolation of figure 7 in 
IKozma &: FranssonI (|1998[ ) to the velocities observed 
in SN 2004dj, we deduced that the hydrogen at the 
observed velocity limits stayed at a temperature of 
T -5000 K up to 859 d and then declined to ~2000 K 
by 1393 d. The f ractional ionization, Xe, was obtained 
from figure 9 of IKozma fc Franss"onl by interpolation. 
This indicated Xe ~ 0.02 at 251 d, falling to Xe ~ 0.0005 
by 1393 d. At each epoch, the T and Xe values were 
assumed to apply to all the hyd rogen within the velocity 
limits. The IKozma fc FranssonI study did not extend to 
epochs as early as the earliest MIR observations (89, 
106, 129 d) of SN 2004dj. We therefore adopted 6000 K 
as a plausible temperature for these epochs and allowed 
the ionization to take the largest value consistent 
with the overall match to the observed continuum. 
Plausible values were obtained, viz: Xe = 0.003 at 89 d 
rising to Xe = 0.005 at 129 d . We used Gaunt factors 
tabulated by [Hummeil (|1988D . For the fb radiation 
we used the co ntinuum recombination coefficients of 
lErcolano fc StOTevI (|2006h . 

The hydrogen maximum velocity, temperature and 
free electron density for each epoch are listed in Table [51 
The fb-ff luminosities are listed in Table [9l col. 4. It can 
be seen that by 500 d the fb-ff contribution to the total 
luminosity is small, becoming increasingly negligible 
at later epochs The fb-ff components for 89-281 d are 
plotted in FigHH At 500 d the fb-ff flux is too weak 
to appear on the 500 d plot which has been scaled to 
allow easy comparison with the earlier epochs. A more 
comprehensive estimation of the fb-ff flux, taking into 



account the temperature and ionization gradients, is 
beyond the scope of this paper. 

As noted in §3.4.3.1, in the wavelength region 
(A > 15 fim) where the cold component makes a 
significant contribution (> 20%) to the total fiux, the 
IS IR echo maintained a near-constant luminosity (15- 
150 fim) of ~ 2.0 X 10^8 ergs-^between 89 d and 1393 d. 
Moreover, in setting the IS IR model to reproduce the 
longwave excess at the latest epochs, when the fb-ff flux 
was negligible, it was found that a satisfactory match 
to the longwave (A > 15 /im) excess was automatically 
achieved for all epochs. Consequently, the fb-ff flux 
could make, at most, only a minor contribution even 
at the earliest epochs. We conclude that any free-free 
emission was too weak to account for the cold component 
at any epoch. 



3.4.4. The Warm Component 
3.4.4.1. Late-Time Optical and NIR Line Profiles 

In §2.4 we described the evolution of optical and NIR 
line proflles in a number of species. In particular, we 
noted the development of a sharp suppression in the 
red wing, suggesting dust formation causing obscuration 
of the far side of the ejecta. Could this dust also be 
responsible for the warm component at later epochs? In 
this and the next sections we explore this possibility. 
Here we examine the distributions of dust that might 
give rise to the late-time optical/NIR line proflles. In 
§3.4.4.2 we shall then test the hypothesis that the 
same dust distributions are responsible for the warm 
component. We modelled the line proflles of Ha, Pa/3, 
[O I] 6300 A, [Fe II] 7155 A and [Fe II] 12567 A. 
The period 461 — 925 d was studied since this was 
mostly overlapped by the 500 — 1393 d during which 
the blackbody analysis and more detailed studies of 
the MIR emission (§3.4.4.2) suggest that substantial 
quantities of dust were present in the ejecta. We note 
(Table |6|) that the Ha line widths towards the end of the 
bright plateau phase greatly exceeded those observed 
on 461 d and later, implying that there was a negligible 
contribution of any light echo of the early phase proflles 
to the late-time proflles. 

We did also examine the blue asymmetry of the Ha 
and [Fe II] 7155 A proflles at a much earlier epoch 
(283 d) but were unable to achieve a satisfactory match 
using the model adopted for later epochs. In any case, as 
already demonstrated in §3.4.2, we were able to account 
for the 251/81 d IR excess as th ermal emission from 
dus t formed in the CDS. Moreover. IChugai et al] ()2005[ ) 
and iChugail (|2006D successfully explained the Ha profile 
up to ~300 d by invoking the emergence of an intrinsic 
asymmetric, bipolar core (i.e. no dust involved). We 
conclude that there is no evidence for the formation of 
new ejecta dust earlier than 461 d. 

The line-profile red wings exhibit increasingly abrupt 
declines after ^1 year (Fig. |TU|) suggesting the conden- 
sation of attenuating dust. In contrast, the extended 
blue wings exhibit little sign of developing suppression. 
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This behavior points to the attenuating dust lying at 
low line-of-sight velocities, such as a face-on disk-like 
distribution centered on the SN center of mass. Nev- 
ertheless, in our initial considerations of possible dust 
configurations we included the case of a spherically 
symmetric dust sphere. Our line-profile model comprises 
a homologously expanding sphere of gas responsible for 
most of the observed line flux, with the emission being 
attenuated by a dust zone lying concentrically with the 
gas sphere. The dust zone is also assumed to participate 
in the expansion. The gas emissivity is assumed to have 
a power-law dependence with radius. 

Several arrangements of dust were examined. These 
were (i) a uniform opaque sphere fFig. 1161 upper panels), 
(ii) a face-on opaque ring, and (iii) a thin face-on disk 
(Fig. 1161 lower panels) whose opacity was uniform or 
varied radially as an power law along the disk plane 
but was always uniform normal to the disk plane. (For 
the disks, "opacity" refers to the value normal to the 
disk plane.) For the radially-varying version of (iii), a 
power law index of 5 = 1.9 was chosen as typical of 
the power law indices determined for the gas emission 
profiles. Configurations (i)-(iii) were initially tested 
against the 895 d optical profiles since these were (a) of 
the highest available resolution, and (b) corresponded 
to a time by which the SN spectral continuum was 
overwhelmingly dominated by thermal emission from 
dust. Configurations (i)-(ii) could not reproduce the 
observed line profiles. This is illustrated in Fig. [16] 
(upper panels) where we show two examples of the line 
profile produced by an opaque, concentric dust sphere, 
compared with the observed [O I] 6300 A profile at 
895 d. Only configuration (iii) - the dust disk - was able 
to reproduce the observed profiles. Within this config- 
uration we examined cases where the dust distribution 
(a) extended beyond the gas limit, (b) ceased abruptly 
at a given velocity within the gas sphere, or (c) extended 
uniformly to a given velocity, Vduni, and then declined 
as a power law to the edge of the gas sphere. Case (a) 
failed to reproduce the observed line profiles. Cases (b) 
and (c) are discussed below. 

As described in §2.4, by 895 d and 925 d, a "sec- 
ondary" , weaker, redshifted peak had also appeared 
in the Ha and [O I] 6300 A profiles, suggesting that 
a fraction of the line flux originated in an emission 
zone centered on the SN and having the geometry of 
an expanding ring, jet or cone. The second peak had 
a redshift of only -1-160 to -1-170 kms~^. Given the 
preference for a thin face-on disk of attenuating dust to 
account for the asymmetry of the line profiles, a natural 
explanation for the low velocity of the second peak is 
that its source was actually at a velocity comparable to 
that of the gas sphere, but moving in a way which was 
roughly coplanar with the dust disk. Consequently the 
low inclination angle led to the low observed redshift in 
the second peak. The simplest geometry which accounts 
for the second peak is a ring-like emission zone whose 
plane lay at a small angle to the "face-on" plane. A cone 
or jet geometry is more problematic. If we insisted on 
maintaining axial symmetry then a cone/jet would have 
to be normal to the disk but with a mysteriously low 



velocity. Alternatively, a cone/jet lying near the disk 
plane, while possibly allowing a high intrinsic velocity, 
would have the unattractive feature of breaking the 
axial symmetry. Therefore, to account for the second 
peak, we added a thin, near-copl anar ring of ernitting 
gas to the hne-profi le model (cf. IGerardv et al.l 120001 : 
iFransson et "al|[2005l . for the case of SN 1998S). The ring 
was assumed to participate in the overall homologous 
expansion of the gas. 

We set the intrinsic ring velocity equal to the fastest 
moving gas observed viz. hydrogen at 2400 kms~^as in- 
dicated by the Ha blue wing on 895 and 925 d, although 
higher velocities could have been used. Consequently 
the ring emission is unattenuated by the dust. For ease 
of computation, we retained a face-on disk (i = 0°). 
Tilting the disk to become exactly coplanar with the 
ring would have had only a small effect on the line profile. 

In matching the model to the observed profiles the free 
parameters for the gas are the maximum (outer) veloc- 
ity of the gas sphere, Vgmax, the gas emissivity scaling 
factor and power-law index fi, the inclination, i^, of the 
ring component assuming an intrinsic expansion veloc- 
ity of 2400 kms~^, and the ring component emissivity 
assumed uniform. In addition, the overall wavelength 
positions of the observed profiles were allowed to vary 
by small amounts (see below). For ease of comparison, 
the dust-disk dimensions were also expressed as velocities 
within the homologous expansion. Thus, for the case (b) 
dust disk the free parameters are the maximum radial 
velocity of the dust disk, Udmax, the maximum velocity, 
Vdthi perpendicular to the disk plane and the magnitude 
of the dust radial density power-law index, 5. Multiplica- 
tion of Wdmax and 2i;dth by the epoch yields, respectively, 
the disk radius and thickness at that time. For case (c) 
i^dmax is replaced with Uduni- We adopted amorphous 
carbon as the grain material (see §3.4.4.2) and a power 
law g rain-size distribution, index m = 3.5 (jMathis et al.l 
Il977f l. with a(jriin) = 0.005 /im and afj^ax) = 0.05 /Ltm. 
The optical depth normally through the disk was then 
calculated. For example, for case (b) the optical depth 
for emission from gas lying behind the disk at wavelength 
A, epoch t and radial velocity Vd, T\{t, Vd), is given by: 

.(3) 

where k is proportional to the grain number density, 
at a fiducial epoch, to the radial velocity and to the 
grain radius, p is the grain material density and k\ is 
the mass absorption coefhcient of the grain material at 
wavelength A. For emission from gas lying within the 
disk the same formula is used but replacing 2wdth with 
the velocity equivalent to the path from the emission 
point to the near side of the disk. The disk thickness was 
always < 20% of the disk diameter and so edge-effects 
at the disk's outer limit were assumed to be negligible. 

The gas emission profiles (i.e., those for the attenu- 
ated sphere and unattenuated ring) were individually 
convolved with the appropriate slit PSF. The PSF 
FWHM values were 6.5, 4.2, 6.5, 1.25, and 1.75 A for 
461, 467, 554, 895, and 925 d, respectively. The two 
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components were then summed to form the final model 
profile which was compared with the observed profile. 
The profile model (case (b) disk) was adjusted to match 
the observed profiles as follows. The parameters I'gmax, 
P and the gas sphere emissivity scaling factor were 
adjusted to provide a match to the blue wing of the 
profile. Then, for a given value of 5 (0 or 1.9), Udth 
(equivalent to the disk half-thickness) together with 
the grain number density scaler, k, were adjusted to 
reproduce the sharp-decline section of the profile. The 
velocity Wdmax was varied to reproduce the suppressed 
red wing. The gas ring emissivity and inclination were 
then adjusted to reproduce the secondary peak. Some 
iteration of all the model parameters was necessary to 
reach the final profile match. For the profile model 
(case (c) disk) Wdmax was replaced with ?;duni, and 5 was 
allowed to vary to arbitrarily large steepness. 

The offset of the sharp-decline section of the profiles 
from zero velocity was only about —100 kms~^ after cor- 
rection for the -1-221 kms~^ of the heliocentric velocity. 
Consequently, the match to the sharp-decline section of 
the profile (i.e., determination of the disk half-thickness) 
could be significantly affected by errors in individual 
wavelength measurements or intrinsic variations in the 
distributions of different gas species. Forcing a match in 
each case by letting the wavelength position of each ob- 
served profile to vary by small amounts, we were able to 
estimate the disk half-thickness and its uncertainty. For 
a heliocentric velocity of -1-221 kms~^ we obtained a disk 
half-thickness equivalent to Wdth = +112 ± 18 kms^^. 
There was no significant variation in Wdth during the 
period covered by the profile study. Note that a change 
in the heliocentric velocity estimate would yield the 
same change in the disk half-thickness velocity. 

As indicated above, we first carried out model matches 
to the 895 d profiles. It was found that the best case (b) 
matches were obtained when the disk was highly opaque 
(say, r > 5) at all radial locations. For case (c), 
comparably good matches were achieved by setting the 
uniform zone at a high opacity (t > 5) and 5 at values 
steeper than ~ 10 — i.e., the dust density beyond the 
uniform zone had to be extremely steep. Less steep 
declines tended to suppress the visibility of the central 
minimum. We conclude that the dust disk was highly 
opaque in the optical region but terminated abruptly 
at an approximately fixed radius (see §3.4.4.2). Given 
that this condition was indicated by both cases (b) 
and (c), we abandoned the more complicated case (c) 
model and completed the analysis for all the 461-925 d 
spectra using only case (b) — i.e., a dust density which 
was uniform or declined radially as r~^'® and which 
terminated abruptly at Wdmax- 

The models for 895 d were initially adjusted to deter- 
mine the minimum opacity - that is, the minimum dust 
mass that was needed to provide a satisfactory match 
to the observed profile. The dust masses were obtained 
from the following: 

8 1 1 

(4) 



At 895 d, with a uniform disk [5 — 0), the 
dust extended to t^dmax = 590 ± 25 kms~^ i.e. 
r-dmax = (4.6 ± 0.2) X 10^^ cm. The disk half-thickness 
was 19% of this. Simultaneous matches to all three 
profiles (Ha, [O I] 6300 A, [Fe II] 7155 A) required 
a minimum r of 13 at 6300 A, corresponding to a 
minimum dust mass of 0.27 x 10^"' Mq. Note that the 
high optical depth was demanded by the very sharp 
decline, unresolved even at the 1.25 A resolution of 
895 d. Other parameters are: r-gmax = 7.65 x 10^^ cm, 
13 = 1.9, (ring flux)/(total flux) = 0.07 and ring tilt 6°. 

In using the dust-disk model to reproduce the contem- 
porary MIR continua (§3.4.4.2) it was found that the 
opacities and dust masses for all lines and epochs had to 
be somewhat higher than the minimum values required 
to match the observed optical and NIR line profiles. 
For example, for [O I] 6300 A at 895 d and a uniform 
disk it was necessary to set r > 18 ± 3, corresponding 
to a minimum dust mass of (0.38 ± 0.06) x lO^'* M0. 
In the rest of this subsection, therefore, we present 
and consider results for profile matching incorporating 
dust-disk masses obtained by interpolating to the profile 
epochs the values obtained from the MIR continuum 
modelling. A sketch of the case (b) [O I] 6300 A 
uniform [5 ~ 0) disk model at 895 d, is shown in 
Fig. [16] (lower LH panel) together with an illustration 
(lower RH panel) of the model match (solid red line) 
to the observed profile (blue). Also shown are the 
intrinsic line-profile contributions from the attenuated 
gas sphere (green) and unattenuated gas ring (cyan) 
components, as well as the final profile (dotted red 
line) which would result if the attenuating dust were 
removed. In Fig. [17] we show all the individual matches 
to the line profiles for 461-925 d. In general, good 
matches were achieved for epochs 554 d, 895 d, and 
925 d. Poorer matches were obtained at 461 d and 467 d. 

On 895 d and 925 d the model did not reproduce the 
extended red wing seen in the [Fe II] 7155 A profile. 
This suggests an additional, faster moving component 
of iron, and may be indicative of ^^Ni asymmetry or 
"bull ets" in the initial explosion (e.g. iBurrows et al.l 
I1995D . At the earliest two epochs (461 d and 467 d) 
the model matches are also poorer. Specifically, if the 
model was adjusted to match the full suppression of 
the extreme red wing, then it also over-attenuated the 
less redshifted portion of the red wing; in other words, 
the steep red decline in the observed profile is generally 
less pronounced at these earliest epochs. Alternatively, 
matching to the less redshifted portion of the red wing 
meant that the full suppression of the extreme red 
wing was not reproduced. These points suggest that, 
during the 461-467 d period, dust formation was less 
complete and was not yet fully opaque over the whole 
disk. We therefore adjusted the models to reproduce the 
suppression of the extreme red wing, but recognise that 
the derived dust mass lower limits for 461-467 d may 
well be overestimated. 

Between 554 d and 925 d the minimum dust mass 
increased from 0.25 x 10"^ M© to 0.4 x 10""* Mq. 
During this time the radius of the dust disk appeared 
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to actually decrease slightly from 5.3 x 10^^ cm to 
4.6 X 10^^ cm. We explored the possibility that this was 
due to a declining optical depth as the disk expanded 
— i.e., that while the dust disk continued to expand, 
the optically-thick/thin boundary declined in radius. 
We found that the requirement of a steep, sharp cutoff 
at the disk limit (see above) ruled out this expla- 
nation implying that the decrease in the dust radius 
was real. This issue will be considered further in §3.4.4.2. 

The above analysis was repeated with S — 1.9. Similar 
disk radii were obtained, but the minimum dust masses 
were ~ x20 larger than for 5 = 0. This is as one 
would expect. To maintain the line-profile matches 
with d = 1.9 required the maintenance of a high optical 
depth at the disk edge in the optical/NIR wavelength 
range. If we demand that the optical depths be the 
same at rdmax for both 6 = and 6 > then it may be 
shown that MdiS)/MdiO) = 2/(2 - 6), where Mdi6) is 
the dust mass for a given value of 6 and Md{0) is the 
mass for ^ = 0. Thus, for 6 = 1.9, Af<j(5)/Md(0) = 20. 
We note that this assumes that the power law continues 
to the center, which is unlikely to be the case. Indeed 
for (5 = 2 or more, the mass would be infinite. A more 
plausible behavior would be that the density law flattens 
toward the center. For example, if we assume that 
within O.OSrdmax the dust density became uniform, with 
6 = 1.9 outside O.OSrdmax, then the total minimum dust 
mass would only be about a factor of 3 larger than for 
the totally uniform disk case — that is, about 10^"* M0 

The model gas sphere expansion velocity remained at 
around 2400 kms~^ in hydrogen at all epochs. Similar 
velocities were obtained in [Fe II] 12567 A at 461 d 
and 554 d, and in [O I] 6300 A and [Fe II] 7155 A at 
467 d. By 895 d and 925 d the expansion velocities in 
[O I] 6300 A and [Fe II] 7155 A had fallen to less than 
1000 kms^^, perhaps suggesting some stratification by 
element. At 554 d and earlier epochs, the ring flux was 
negligible. At the later epochs, the ring inclination is 
6° (i.e., close to face-on). The fractional contribution 
of the ring to the total observed flux never exceeded 
20%. Indeed, were it not for the dust disk, it is unlikely 
that the relatively weak ring emission would have been 
detected since it would have been swamped by the 
emission from the gas sphere (see Fig. [T5|. We therefore 
regard the ring component as an interesting but minor 
effect. 

We conclude from the line-profile analysis that, in 
the 461-925 d period, the mass of the dust disk and its 
radial extent can be small enough to be consistent with 
an origin as newly- condensed ejecta dust. In §3.4.4.2 we 
shall show that the same dust disk can account for the 
MIR continua. 



3.4.4.2. Later Phase IR Excess and Ejecta Dust Conden- 
sation 

We have argued that, up to at least 281 d, the warm 
component could be fully accounted for in terms of an 
IR echo from a CDS. By 500 d the warm component 



luminosity exceeded the CDS contribution by a factor 
of ^ 20. This "late IR excess" , together with the steady 
shift of the peak emission of the SED to longer wave- 
lengths in the period ~ 250-1200 d (Figs. HIS]) suggests 
the appearance of an additional population of warm, 
but cooling dust. What is the location, distribution 
and heating mechanism of this dust? Such IR emission 
can arise from (i) heating, by a number of possible 
mechanisms (see below), of new dust formed either in 
the ejecta, or in a late surge of dust growth within the 
CDS or (ii) heating of pre-existing circumstellar dust by 
the early-time SN luminosity (a conventional IR echo), 
or possibly the forward shock travelling into undisturbed 
CSM beyond the CDS. We can dismiss immediately 
hypothesis (ii): 

(a) An IR echo does not naturally produce the long 
delay (at least ~ 300 d) before the commencement of 
the warm component excess flux rise. To account for 
this within an IR echo scenario, we would have to invoke 
an ad hoc asymmetry in the CSM. 

(b) IR-echo or forward-shock heating would not 
account for the observed late-time red wing suppression 
in the spectral line profiles discussed in §3.4.4.1. In 
contrast, red-wing attenuation can easily be produced 
by dust formation in the ejecta. New CDS dust could 
conceivably also have produced such an effect provided 
that the dust did not lie completely outside the opti- 
cal/IR emission zone. 

(c) Pre-existing CSM dust, whether heated by the 
SN luminosity or a forward shock, would not produce 
steepening of the optical-NIR decline rate simultane- 
ously with the appearance of the late IR excess. In 
contrast, new dust in the ejecta or CDS can be indicated 
by optical-NIR steepening. In Fig. SI after ~470 d the 
VRIJ light curves show weak evidence of steepening, 
providing a minor additional argument against a CSM 
IR echo origin for the later warm component, although 
the errors on the latest points are large. We also note 
that an alternative explanation for such steepening 
could be a faster-than-expected dec rease in the 7 -ray 
absorption relative to that invoked bv lLi et all ()1993 ^ for 
the radioactive energy deposition in SN 1987A (Fig. |4]). 

The above points, especially (a) and (b), leaves us 
with hypothesis (i) viz. that the late IR excess is due to 
emission from new dust formed in either the ejecta or 
CDS. The later-era appearance of the red wing suppres- 
sion coincides roughly with the emergence of the late IR 
excess, suggesting that the same dust could have been 
responsible for both the optical/NIR attenuation effects 
and the rise of the MIR emission. That this dust formed 
in the ejecta rather than the CDS is indicated by the 
fact that the extent of the attenuating dust derived from 
late-time line-profile analysis was comparable to that 
of the MIR-emitting dust as derived from blackbody 
analysis; for example, at 859 d the blackbody radius is 
3.3 X 10^^ cm compared with a profile-derived radius of 
(4.4 ±0.2) X 10^^ cm at 895 d. Moreover, the line-profile 
analysis demonstrates that the dust was of high opacity 
in the optical/NIR, and distributed as a near face-on 
disk lying concentrically with the SN center of mass. 
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Attenuation by a CDS could not have produced the 
late-time line profile behavior. 

In view of the above points, we shall now give detailed 
consideration to the scenario where the warm compo- 
nent emission originated in newly-formed ejecta dust. 
Such dust could be heated by a number of mechanisms 
including radioactive decay, ionization freeze-out effect, 
embedded pulsar or reverse-shock radiation. In the 
case of SN 2004dj we can rule out radioactivity as the 
principal heating mechanism of any putative new ejecta 
dust. Inspection of Table [7] shows that, as early as 
500 d, the warm component flux exceeded the total 
deposited radioactive luminosity by a factor of 2 growing 
to a factor of ~ 450 by 1393 d. Indeed, by the latest 
epoch, the warm component luminosity was still as high 
as ~ 4 X 10^* ergs~^. We also ru le out the freeze-out 
effect as the main energ y source (jClavton et al.l [l992l : 
iFransson fc Kozinalll993[ ) since, at ~ 1400 d, its contri- 
bution would only be '--^ 5 x 10'^^ ergs. 

Power input from an embedded puls ar, via a pulsar 
wind nebula, is another pos si bility ( Wooslev et al.lll989l : 
iChevalier k, FranssonI IT992I ) . iChevalier fc: FranssonI pre- 
dicted that a distinctive feature of such a source would 
be the presence of certain high ionization lines. The 
earliest epoch they studied was at 1500 d. The nearest 
of our optical/NIR spectra to this epoch are those at 
895 d and 925 d. Of the high ionization lines for which 
IChevalier fc FranssonI make luminosity predictions, only 
[O III] 4959/5007 A was covered by these spectra. We 
co-added the spectra at the two epochs and searched for 
the stronger [O III] 5007 A component. A broad, low 
S/N emission feature was detected with a dereddened 
luminosity of (0.8 ± 0.2) x 10^^ ergs"^ This is about 
half the luminosity predicted by IChevalier fc FranssonI 
A serious difficulty is that the redshift-corrected position 
of the feature lies about -1-400 kms~^ from where it 
would be if due to an [O III] 5007 A line subject to the 
same dust attenuation deduced in other lines (§3.4.4.1). 
Even without dust attenuation, the feature would still be 
-1-200 kms~^ too far to the red. This, together with the 
low S/N, leads us to conclude that there is no persuasive 
evidence for the presence of pulsar-driven high-ionization 
features in the latest spectra of SN 2004dj. This leaves 
us with reverse-shock heating of the ejecta following the 
ejecta-CSM collision as the most promising mechanism 
for the late-time energy source 

To test further the hypothesis that the same dust 
was responsible for the late-time line-profile red-wing 
suppression and the MIR emission we modelled the 
IR emission over a range of late-time epochs using the 
same dust-disk configuration as was employed in the 
line-profile analysis. The emission was derived following 
a similar proce dure to that used in th e isothermal 
dust model by iMeikle et all H20Q7f ) and iKotak et all 
(1^091). The resulting flux was then added to the other 
continuum components and the net model compared 
with the observed continua. 

We considered the thermal radiation from a warm, 
isothermal, face-on disk of dust located symmetrically 



about the SN center of mass. We shall refer to this as 
the IDDM (isothermal dust-disk model). The disk radius 
is f'dmax and, as in §3.4.4.1, the dust number density de- 
clines as r'^^ where 5 is set as or 1.9. As mentioned 
above the disk thickness was < 20% of the disk diame- 
ter and so we judged that a thin-disk treatment would 
provide an adequate means of estimating the flux from 
the disk (i.e., edge effects are ignored). The observed 
IDDM fiux, dF\{t,rd), at wavelength A and time t from 
an elemental ring lying between radii and + dr^ is: 

dFx{t,rd) = 2nrddrdD-^Bx{t){l-exp{-Tx{t,rd))) (5) 

where D is the distance of the SN, B\{t) is the Planck 
function, and T\(t,rd) is the optical depth perpendicu- 
larly through the disk. The total fiux is then found by 
integrating from = to = rdmax- The dust mass is 
obtained from equation (4). Amorphous carbon grains 
were assumed (silicate grains are discussed below). As 
explained before, during the 89 — 281 d period ejecta 
dust formation was unlikely and in any case the early 
IR excess was explainable as emission from CDS dust. 
Therefore the IDDM was introduced at 500 d and used 
at all subsequent epochs. By 652 d the CDS IR emission 
was negligible and so was not included in the modelling 
of this or subsequent phases. We first describe the 
results with J = in the IDDM. The matches together 
with the observations are shown in Figs. [TBI and 1191 The 
overall model parameters are listed in Table IH 

The three MIR observation epochs 500 d, 652 d and 
859 d lay within the timespan of the line-profile analysis. 
We therefore imposed the constraint that the disk radius 
had to be consistent with those values derived from 
the line profiles - that is, rdmax ~ 5 x 10^^ cm and 
'''dth = "^dth X t where Vdth = -1-112 kms~^, with k (and 
therefore Md) being constrained by the demand that 
the disk have a high optical depth in the optical/NIR 
region (§3.4.4.1). In practice, matching the IDDM to 
the MIR continuum demanded optical depths which 
were higher than the minimum values obtained from 
the line-profile analysis. The disk parameters for the 
specific MIR epochs were set by linear interpolation of 
the line-profile-derived values for Tdmax and rdth- The 
k parameter could take values at or above those set 
in the line-profile analysis. Only Td was a completely 
free IDDM parameter. By 500 d the MIR continuum 
up to ^^20 ^m was overwhelmingly due to the warm 
dust disk, with the CDS component yielding no more 
than ~8% of the total flux at any wavelength. For 
A > 20 /im, the IS IR echo dominated. Also from 500 d 
onwards the hot continuum and fb-ff contributions to 
the MIR were negligible. In other words, from 500 d 
onwards, only the dusty disk and IS IR echo made 
significant contributions to the total MIR continuum. 
Consequently, and in order to show the MIR behavior 
in more detail, the optical/NIR region is not shown in 
the plots for 652-1393 d (Fig.[l9l). 

In spite of the line-profile constraints on rdmax and 
Tdth, for epochs 500 d, 652 d, and 859 d we were 
nevertheless able to obtain fair matches to the observed 
continuum (Figs. [T8l 119)) . It was found that the matches 
required the dust to be optically thick (r > 1) in the 
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optical/NIR region, but optically thin in the MIR 
region. A blackbody spectrum matched to the shorter 
MIR wavelengths overproduced the flux at longer 
wavelengths. This is illustrated in Fig. [18] (500 d) and 
Fig. [19] as dotted cyan lines. This restriction, together 
with the fixed values for rdmax and rath allowed us to 
obtain specific values, not limits, for T^, r and Md- At 
500 d the ejecta dust mass was (0.22 ±0.02) x lO""* Mq, 
with Td = 650 ± 15 K. The optical depths through 
the disk were ~10 in the V band and 0.45 ± 0.05 at 
10 /im. By 859 d the dust mass had increased to 
(0.33 ± 0.05) X 10-" Mq, while Td feh to 570 ± 15 K. 
A particular value of this 500-859 d study is that it 
demonstrates that the same dust-disk parameters can 
account for the line profiles and the MIR continuum. 
This adds considerable weight to our contention that 
the two disks are one and the same. Moreover specific 
values, not limits, for the dust masses and temperatures 
were determined for the 500-859 d period. 

As the SN aged between 500 d and 996 d it was found 
that the IDDM steadily approached the blackbody case. 
Moreover, there were no line profiles available after 
925 d. Linear extrapolation to 996 d of the line-profile- 
derived values for rdmax gave rdmax ~ (4 ± 1) x 10^^ cm. 
The model match yielded dust with a high optical depth 
in the optical/NIR but with r - 1 in the MIR. The 
dust mass was Md - 0.5 x lO"'' M© and Td = 520±30 K. 

By epochs 1207 d and 1393 d we found that the 
MIR continuum was best reproduced by allowing a 
continuing increase in the IDDM optical depth such 
that the disk was optically thick at all observed wave- 
lengths (see Table [S]). Indeed, given the uncertainties 
in the observed fluxes, we cannot rule out a totally 
opaque disk at all wavelengths covered. Consequently, 
the IDDM could provide only dust mass lower limits 
of 1.0 X 10""* Mq at 1207 d and 1.5 x lO"'^ Mq at 
1393 d. The 1207 d and 1393 d plots shown in Fig. [H 
are with rdmax and Td set at the limiting values (Table [5]). 

In Table [9] we show the luminosities of the IDDM {S — 
0) components compared with the radioactive deposition 
power of 0.0095 Mq of ^^Ni, including ^^Ni decay. The 
IS IR echo component is excluded since it is powered 
primarily by the SN peak luminosity. The luminosity of 
the thin disk, L/_d_dm, is approximated by: 

LiDDM ~ 27ri?^7rB,(l - exp{-2T,)). (6) 

As already indicated by the BLC analysis (§3.2) the 
post-30 d evolution of SN 2004dj falls into three phases. 
At 89 d, (about half-way down the plateau-edge) the 
luminosity is still dominated by the shock-ionized ejecta, 
with radioactive decay contributing a small proportion 
of the total. There is then a short period just after the 
end of the plateau when radioactive decay deposition 
dominated the luminosity. This is supported by the fact 
that at 129 d, the radioactive luminosity, Lrad, is highly 
similar to the total continuum model luminosity, Ltotai- 
At the nebular phases of 251 d and 281 d, the small 
excess in Lrad, relative to Ltotai, presumably went into 
powering the line emission which is not included in the 
model. By 500 d, in spite of the exclusion of much of 



the line luminosity from the model, Ltotai is x2.5Lrad-, 
rising to x200Lrad by 1393 d. As argued above, the 
most plausible additional power source available at this 
stage is a reverse shock. 

Model matches were also carried out with (5 = 1.9 in 
the IDDM component. Similar dust parameters were 
found to those obtained for ^ = 0. The only significant 
difference was that the dust masses were x3 larger, 
assuming a uniform density distribution within the inner 
0.05rdmax (cf. §3.4.4.1). These larger values are due to 
the growing proportion of dust mass concentrated in the 
optically-thick region of the disk. 

We also considered a dusty disk of warm silicate 
grains, with the radii and thickness determined by the 
line-profile analysis as before. The demand that the disk 
should be of optical depth in the optical/NIR up to its 
edge meant that the replacement of amorphous carbon 
dust with silicate dust had no effect on these dimensions. 
The problem with silicate dust continuum matching 
is the absence of the 8 — 14 /im silicate feature in the 
observed continua during the period when the dust was 
optically thin in the MIR, up to 996 d. Attempts to 
suppress the feature in the IDDM by increasing the dust 
mass and hence the optical depth yielded a continuum 
that was too bright. Only for 1207 d and 1393 d was 
silicate dust able to reproduce the observed continuum. 
This is not surprising since, by these two epochs, the 
warm ejecta dust was close to being opaque over the 
MIR range observed. Between 500 d and 996 d the 
proportion of silicate grains by mass could have been 
no more than 20%, and was usually significantly less 
than this. Indeed, the data were always consistent 
with there being no silicate grains at all. We also note 
the absence of the SiO feature at 7.5 — 9.3 ^m. SiO 
formati on is a necessary step on the way to silicate 
grains (jTodini fc Ferrarall200ll: iNozawa et al.ll2003D . 

The points made in the previous paragraph argue 
against a substantial amount of silicate dust in the ejecta 
of SN 2004dj. The absence of silicate grains is consistent 
with the presence of strong CO fundamental and first 
overtone emission in the period ~100 d to 300 — 500 d. A 
high C/0 ratio in the ejecta could result in most of the 
oxygen being absorbed as CO leaving behind an excess 
of carbon to provide carbon grains, but little oxygen to 
provide silicate grains. However, in the SN environment, 
the net grain population can be affected by factors 
in addition to the C/O ratio. These include molecule 
destruction by high-energy electrons, charge transfer 
reactions and ejecta density (.L iu & Dalgarno 199^ 
INozawa et all 120031 : iDeneault et al.ll2006l ). Nevertheless, 
we conclude that the dust grains in SN 2004dj were 
predominantly composed of non-silicate material. 

The success in reproducing the MIR continua using 
the same dust disk as was invoked to explain con- 
temporary line profiles tends to support the curious 
result from §3.4.4.1 that, rather than expanding, the 
radius of the dust disk actually shrunk by a small 
amount. Indeed, if we include the disk radii for 1207 d 
and 1393 d derived from the IDDM matches (Table [H] 
col. 7) we find a shrinkage of 27% since 500 d. In 
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SN 2004et ()Kotak et alj I2OO90 it was found that the 
dust radius remained roughly constant. The explanation 
offered was that the dust was contained within an 
optically-thin cloud of optically-thick, pressure-confined 
clumps. But the availability of late-time optical/NIR 
spectra for SN 2004dj and the profile modelling pre- 
sented here shows that a clumping explanation for the 
non-expanding dust-disk radius would not work. We 
note that, during the 500-859 d period in SN 2004dj, 
the product of the blue-wing, half-maximum (BHM) 
velocities and epoch for the line profiles, Rbhm, yields 
a roughly constant value of ~ 4 x 10^^ cm (see Table [6]) 
- similar to that obtained for the dust-disk radius. This 
coincidence may imply that the extent of both the dust 
and the bulk of the ejecta gas emission was physically 
constrained within a radius of ~ 5 x 10^^ cm. 

We suggest that the apparent shrinkage of the dust- 
disk radius may have been due to the presence and 
inward motion of the reverse shock such that the shock 
position defined the disk radius. Dust formation could 
have continued within the disk with the dust taking part 
in the overall ejecta outfiow, but as the ejecta passed 
through the reverse shock it would have been destroyed. 
There would have been a net increase in dust mass with 
time, as observed, if the dust growth rate within the 
disk exceeded the destruction rate. A possible difficulty 
with this scenario is that the dust mass appeared to 
continue to grow right up to the final epoch at 1393 d. 
This is rather later thaii dust formation studies suggest. 
iTodini fc FeTraral (|200l[ ) find that all grain condensa- 
tion would be co i nplete by 8 00 days. Moreover both 
ITodini fc FeTraral (|200l and iNozawa et all (|200l find 
carbon dust condensation is complete in not much more 
than one year. An alternative explanation, therefore, 
might be that the dust mass in SN 2004dj did not in- 
crease after ~400 d, but rather that the density gradient 
of the outer region of the disk was actually steeper than 
r~^. Thus, as the ejecta expanded, the dust density at 
a fixed location in the outer region would have grown, 
yielding an apparently higher total dust mass — i.e., 
more dust emerged from the optically-thick inner regions. 

Radiation from the reverse shock could have been 
primarily responsible for heating the dust. As more 
and more ejecta passed through, the shock and its 
radiation would have weakened causing the dust to cool, 
as observed. The reverse shock may also have been 
responsible for the approximate constancy of the radius 
of the line emitting gas sphere. Further examination 
of the reverse-shocked dust-disk hypothesis is beyond 
the scope of this paper. We note th at in a recent 
observ ational study of Cassiopeia A iDelanev et al.l 
(|2010( ) deduce a fiattened ejecta distribution or "thick 
disk" containing all the ejecta structures. They also 
deduce the occurrence of a roughly spherical reverse 
shock. 

Intrinsic axial asymmetry in the form of a bipolar 
jet biased towards the observer has been invoked by 
iChugai et al.l (|2005[ ) to account for the Ha profile in 
SN 2004dj up to about 1 year. We note that the jet 
makes an angle of only ^ 15° to the normal to the dust 
disk plane derived from our line-profile analysis, perhaps 



indicating a physical connection. However, line-profile 
asymmetry in the earlier (pre-'^l year) nebular spectra 
of other CCSNe also tend to be blue-biased. This implies 
that, in general, an intrinsic bias towards the observer 
cann ot be the explanation for such line blueshifts (cf. 
iMilisavlievic et al.. 2010. 1. At later (post-'-^^l year) nebu- 
lar ep o chs, line blueshifts arealso often seen (iLucv et al.l 
IQSgt iSpvromilio et al.l 119901: iTuratto et al.l 119931 
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20001: iFassia et ai.l 120021: lElmhamdi et al.l 120031: 

Pozzo et al. I I2004n and these are usually attributed 
to attenuation by newly-formed dust in the ejecta or 
CDS. Our dust-disk model invokes this scenario. This 
has the advantage of explaining the red-wing sup- 
pression without invoking an intrinsic observer-biased 
axial asymmetry in the SN. Moreover, our model uses 
the same dust disk to simultaneously account for the 
line-profile attenuation and t he MIR emi ssion. We 
consider it unlikely that the IChugai et al.l model can 
provide a superior alternative explanation for the line 
profiles presented and analysed in §3.4.4.1. Further 
examinati on of the r elationship between our model and 
that of Ch ugai et al.l is beyond the scope of this paper. 

We conclude that dust formation in the ejecta of 
SN 2004dj had commenced, and may even have been 
completed, by 500 d with a near-face-on, disk-like 
distribution. This dust was responsible both for the 
late-time line-profile red-wing suppression and the bulk 
of the MIR luminosity up to ^ 20/im. The main source 
of dust heating was probably reverse-shock radiation. 
Assuming 5 = 0, the dust mass was at least 10"'' Mq. 
For 6 = 1.9, the lower limit rises to about 3 x 10~* Mq. 
The value of 6 is poorly constrained. We reject silicates 
as the grain material. 

4. CONCLUSIONS 

We have presented optical, NIR, and MIR observa- 
tions of the Type IIP SN 2004dj. The combination of 
wavelength and temporal coverage achieved makes this 
SN one of the most closely studied of such events. In the 
present work we have analyzed the SN continuum over a 
period spanning 89-1393 d, augmented by a line-profile 
analysis over 461-925 d. Our conclusions are as follows. 

(1) A mass of 0.0095 ± 0.002 Mq of ^^Ni was ejected 
from SN 2004dj, which is less than that reported by 
other authors. The period during which the radioactive 
tail dominated the bolometric light curve lasted for an 
unusually short period of only ~35 d. Subsequently, a 
different energy source dominated; we suggest reverse- 
shock heating. 

(2) At early times the optical/NIR ("hot") part of 
the continuum provided most of the SN luminosity. 
This emission is attributed to hot, optically-thick ejecta 
gas. At later epochs the optical/NIR luminosity was 
increasingly due to nebular emission, and formed only a 
minor, ultimately negligible, proportion of the BLC. 

(3) At both early and late times, the long-wave 
portion of the MIR continuum (the "cold" component) 
was primarily due to an IS IR echo. Free-free radiation 
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m ade only a minor contr ibution. As originally suggested 
by iBode fc EvansI (|1980D . the analysis of SN IR echoes 
may provide a useful way of studying IS dust in nearby 
galaxies. 

Subse quent to the submission of this paper, iSzalai et al.l 
(|2011l ) reported their findings on SN 2004dj, making use 
of some of the data presented in this work. They dismiss 
the IS IR echo, and argue against pre-existing dust since 
the estimated extinction of SN 2004dj was lower than 
that of other SNe. We do not concur. We have shown 
that an extinction of only Av=0.02 i s all that is required 
to account for an adequate IR echo. ISzalai et al.l (|2011[ ) 
also argue that the early UV/X-ray flash would create 
a dust-free cavity of up to 10^^ cm and that OB-stars 
would have expelled the ISM/dust from the cluster. In 
fact most of the IS IR echo longward of ~ 10 /im comes 
from dust lying at considerably more than 1 parsec. 
We therefore strongly favour a scenario whereby the 
long- wave component was primarily due to an IS IR echo. 

(4) The early-time NIR/MIR ("warm") component 
was probably due to thermal emission from non-silicate 
dust formed in the CDS. The CDS dust growth began 
at about 50 d, reached 90% of maximum by 165 d, 
and approached a maximum of 0.33 x 10~^ M©; the 
dust mass produced in this way was small. Heating 
of the dust by the contemporary optical-NIR BLC 
completely accounts for the streng th and evolu t ion o f 
the early-time warm component. ISzalai et al.l (|2011 f) 
did not present contemporary optical or NIR data and 
consequently did not identify the early warm component. 

(5) The late-time warm component was dominated 
by the luminosity of newly-formed, non-silicate dust in 
the ejecta. The dust growth commenced at some time 
between 281 d and 461 d. The same dust was responsible 
for the late-time red-wing attenuation of optical and 
NIR spectral line profiles. The dust was distributed as a 
nearly face-on disk, within a spherical cloud of emitting 
gas. During 500-996 d the disk was effectively opaque 
in the optical/NIR region, but was optically thin at 
longer wavelengths. The dust mass appeared to grow 
during this period, attaining (0.5 ± 0.1) x 10~* Mq by 
996 d for a uniform density disk, or a few times more 
than this for an r~^'^ gradient outside 0.05rdmax and 
uniform within. However, it may be that the dust mass 
"growth" was really due to the emergence of previously 
formed (i.e., pre-500 d) dust from optically-thick regions 
of a disk with an even steeper density gradient. For 
the latest two epochs (1207 d and 1393 d) the dust 
was optically thick at all wavelengths and only lower 
limits could be obtained for a given gradient — for 
example, > 10~^ Mq for a flat gradient and a factor of 
3 higher limit for an r~^-^ gradient outside 0.05rdmax 
and uniform within. This is broadly in agreement with 
ISzalai et al.l (pOll.) . For a smooth distribution, they 
found a dust mass of ~ 10~^ — 10~* Mq though, with 
clumping, up to ~ 10~^ Mq would also be possible. 

(6) Rather than expanding, the dust-disk radius 
appeared to slowly shrink. This may have been due to 
the dust extent being confined by the reverse shock. 



which also heated the grains radiatively. 

(7) While the latest epochs provide only lower limits 
to the mass of dust produced by SN 2004dj, these limits 
are at least a factor of 100 below the 0.1 Mq of grains 
per SN required to account for the dust observed at 
high rcdshifts. Moreover, measurements as late as 996 d 
yield actual dust masses of only - 10-* Mq. While not 
completely ruling out the possibility that typical CCSN 
ejecta are major contributors to cosmic dust production, 
this study does suggest that, at least for SN 2004dj, the 
dust-mass production was small. 
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TABLE 1 
MiD-IR Photometry of SN 2004dj. 



Flux (mjy)** 







Epoeh* 


tcxp 




IRAC 




PUI 


IVilr^o 


Date 


MJD 


(d) 


(s) 


3.6 


4.5 /^m 


5.8 /^m 


8.0 fiiii 


16 )Um 


22 /im 


24 /^m 


zUU4 (Jet U7 


53285.1 


89.1 


134 


11.25(2j 


10.00(2) 


6.85(3) 


4.32(4) 








2004 Oct 08° 


53286.4 


90.4 


107 


10.94(2) 


10.01(2) 


6.74(3) 


4.09(5) 








2004 Oct 12" 


53290.1 


94.1 


107 


8.30(2) 


9.71(2) 


6.02(3) 


3.45(4) 








2004 Oct 14" 


53292.5 


96.5 


160 














0.74(16) 


2004 JNov 01 


53310.1 


114.1 


134 


4.37(2) 


7.70(2) 


4.78(3) 


2.49(3) 








2004 Nov 06" 


53315.7 


119.7 


160 














0.70(9) 


2005 Mar 03" 


53432.3 


236.3 


160 














0.89(9) 


2005 Mar 24" 


53453.9 


257.9 


134 


1.25(2) 


3.50(1) 


1.57(3) 


1.27(4) 








2005 Apr 02" 


53462.0 


266.0 


160 














0.75(17) 


2005 Oct 20^^ 


53663.7 


467.7 


536 


1.86(2) 


2.36(1) 


2. 4:2(2) 


2.36(4) 








2Uuo INov 22 


53693.9 


497.9 


630 










1.23(4) 






2UUD Mar 2o 


53817.9 


621.9 


536 


1.34(2) 


1.74(1) 


2.38(2) 


2.58(4) 








2UUD Apr Uz 


53827.6 


631.6 


10 














1.0(2)T 


2006 Apr 2d 


53848.4 


652.4 


630 










1.43(7) 






2006 Oct 28 


54036.2 


840.2 


536 


0.85(2) 


1.14(1) 


1.81(2) 


f\f\ ( A\ 

2.09(4) 








2006 Oct Sl-'^ 


54039.1 


843.1 


14 


0.83(2) 


1.13(2) 


1.50(7) 


2.07(4) 








zUUo JNov lo-' 


54054.0 


859.0 


57 










1.38(7) 






2006 JNov 16 


54055.0 


859.0 


315 










1.47(4) 


1.15(3) 




2006 Dec Ol-* 


54070.3 


874.3 


160 














1.15(15) 


2007 Mar 24*^ 


54183.0 


987.0 


1258 










1.28(9)'"'' 


0.96(8)1"' 




2007 Mar 24-' 


54183.0 


987.0 


132 










1.37(10)1 T 






2007 Apr 02^ 


54192.8 


996.8 


536 


0.50(2) 


0.67(1) 


1.15(2) 


1.37(4) 








2007 Apr 02-' 


54192.8 


996.8 


125 


f\ A(\ { c\\ 

0.49(2) 


0.67(1 ) 


1.13(3) 


1.38(4) 








OAn'7 A 1 Q f 

200/ Apr ivS-' 


O420o.l 


1007.1 


160 














o.yo^ioj 


2007 (Jet 24^ 


54397.3 


1201.3 


494 














1.09(13) 


2007 Oet 24^ 


54397.3 


1201.3 


1360 














1.07(13) 


2007 Nov 04^^ 


54408.3 


1212.3 


283 










1.18(4) 






2007 JNov 19^ 


54423.3 


1227.3 


322 


0.36(2) 


0.45(1) 


0.76(2) 


1.01(4) 








2007 Nov 23^ 


54427. 1 


1231.1 


322 


0.36(2) 


0.44(1) 


0.76(2) 


1.03(3) 








2007 Nov 24'' 


54428.0 


1232.0 


3485 


0.36(1) 


0.45(1) 


0.80(2) 


1.02(3) 








2007 Nov 29^ 


54433.4 


1237.4 


494 














1.20(6) 


2007 Dec Ob-' 


54440.0 


1244.0 


283 










l.OU^^oO ) ' ' 






zuU( uec lo 


04I:40U.U 




1 ofin 
izDU 










1.12(3) 






2008 Mar 25'' 


54551.0 


1355.0 


1260 










1.03(6) 






2008 Apr 07'' 


54563.8 


1367.8 


3485 


0.33(2) 


0.38(1) 


0.67(3) 


0.90(4) 








2008 Apr 123 


54568.2 


1372.2 


322 


0.33(2) 


0.38(1) 


0.63(2) 


0.95(4) 








2008 Apr 143 


54570.4 


1374.4 


494 














1.24(9) 


2008 Apr 15'' 


54571.0 


1375.0 


1.360 














1.11(14) 


2008 May 039 


54589.3 


1393.3 


283 










1.07(6) 






2009 Nov 14' 


55149.9 


1953.9 


3276 


0.31(1) 














2009 Dec 03' 


55168.7 


1972.7 


3276 


0.30(1) 














2010 Apr 16' 


55302.7 


2106.7 


3276 


0.35(1) 














2010 May 23' 


55339.5 


2113.5 


3276 


0.31(1) 














Saiidagt! 9() 








0.28(."')) 


0.21(1) 


0.11(3) 


0.08(2) 


0.021(7) 


0.013(1) 


0.011(3) 



Note. — * We assume an explosion date of 2004 Jul. 10.0 (MJD=53196.0). 
** All measurements were made in a ~ 3'.'7 radius circular aperture, with annular sky measured between xl.5 and x2.2 the aperture 
radius. The aperture was centered on the WCS co-ordinates of the SN. Statistical uncertainties in the last one or two significant figures 
are shown in brackets. The fluxes have not been corrected for reddening, nor for the contribution of the S96 cluster. The last line shows 
the estimated contribution to the total flux by S96 (see §3.1). The exposure times in Column 4 are per band. It can be seen that some 
same-wavelength observation epochs have negligible temporal spacing as well as differing exposure times (viz. 840/3 d, 859 d, 987 d, 996 d, 
1201 d, 1227/31/32, 1244/54, 1367/72, 1374/5). This was due to the impact of Spitzer scheduling constraints on the MISC and SEEDS 
programs. 

"FID. 00226 Van Dyk (MISC); 
''FID. 00159 Kennicutt et al.(SINGS); 
=FID. 20256 Meikle et al.(MISC); 
'^FID. 20321 Zaritsky; 
<=FID. 30292 Meikle et al.(MISC); 
■''FID. 30494 Sugerman et al.(SEEDS); 
SFID. 40010 Meixner et al. (SEEDS); 
''FID. 40619 Kotak et al.(MISC); 
'FID. 61002 Preedman et al. 

t Target ~50% off edge of field. Total flux estimated by extrapolation. 
Contamination to north-east. 
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TABLE 2 
Near-IR Photometry of SN 2004dj 



Date 


MJD 


Epoch* 


tcxp 




Magnitudes** 




(UT) 




(d) 


(s) 


Z 


J 


H 


K, 


2004 Nov 25 


53334.0 


138.0 


J{70) H{135) 




13.56(4) 


13.24(4) 




2005 Mar 17 


53446.0 


250.0 


Z{60) J(60) H{60) K(50) 


15.4(1)1 


14.35(4) 


14.15(6) 


14.13(5) 


2005 May 11 


53501.2 


305.2 


J(300) H(225) Ks(180) 




14.88(1) 


14.54(2) 


14.42(4) 


2005 May 12 


53502.1 


306.1 


J(300) H(225) 




14.88(1) 


14.52(2) 




2006 Jan 14 


53749.9 


553.9 


J(675) H(1025) Ks(1215) 




15.87(5) 


15.19(5) 


14.75(5) 


Sandago 96 










15.93(11) 


15.70(12) 


15.15(14) 



Note. — The 2005 May 11,12 images are from the TIFKAM IR camera on the 2.4m Hiltner Telescope of the MDM Observatory. All 
the other images are from the LIRIS IR imager/spectrograph on the 4.2m William Herschel Telescope of the Observatorio del Roque de 
los Muchachos, La Palma, 

* We assume an explosion date of 2004 Jul. 10.0 (MJD=53196.0). 

** All measurements were made in a 3.7" radius circular aperture, with annular sky measured between xl.5 and x2.2 the aperture radius. 
The aperture was centered on the centroid of the SN. The error on the last one or two figures is given in parentheses. Errors shown include 
uncertainties in the magnitudes of the four 2MASS comparison field stars. The fluxes have not been corrected for red dening, nor for the 
contribution of the S96 cluster (see §3.1). The last line shows magnitudes of S96 measured from 2MASS images iSkrutskie et al.|[T997f) . 
1 Estimated by extrapolation of JHKs SEDs of calibration field stars to Z band. 



TABLE 3 

MiD-IR Spectroscopy Log of SN 2004dj. 



Date 


MJD 


Epoch* 


texp 


Program 




(UT) 




(d) 


(s) 








2004 Oct 24 


53302.32 


106.3 


610 


PID. 


00226 


Van Dyk (MISC) 


2004 Nov 16 


53325.77 


129.8 


610 


PID. 


00226 


Van Dyk (MISC) 


2005 Mar 18 


53447.65 


251.6 


610 


PID. 


00226 


Van Dyk (MISC) 


2005 Apr 17 


53477.92 


281.9 


610 


PID. 


00226 


Van Dyk (MISC) 


2005 Nov 22 


53696.83 


500.8 


1828 


PID. 


20256 


Meikle et al.(MISC) 


2006 Apr 23 


53848.47 


652.5 


1828 


PID. 


20256 


Meikle et al.(MISC) 


2006 Nov 16 


54055.08 


859.1 


3657 


PID. 


30292 


Meikle et al.(MISC) 


2007 Oct 30 


54403.45 


1207.4 


117 


PID. 


30494 


Sugerman et al. (SEEDS) 


2008 May 03 


54589.35 


1393.3 


3657 


PID. 


40619 


Kotak et al.(MISC) 



Note. — Spectral ranges at each epoch were Short-Low (SL) first and second orders: 7.4 — 14.5 and 5.2 — 8.7 ^m, respectively, 
tcxp gives the exposure time per order, in seconds. 
* We assume an explosion date of 2004 Jul. 10.0 (MJD=53196.0). 
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TABLE 4 

Optical and Near-Infrared Spectroscopy Log of SN 2004dj. 



Date 




MJD 


Epoch* 


Spectral range 


Telescope 


(UT) 






(d) 


(/xm) 




2004 


Oct 05 


53284.9 


88.4 


0.420-0.766 


DDO 1.88m 


2004 


Oct 13 


53291 


95 


0.480-0.910 


Lick 3m 


2004 


Oct 18 


53296 


100 


0.370-0.750 


LAT Scorpio 


2004 


Nov 15 


53324.21 


128.2 


0.402-0.839 


DDO 1.88m 


2004 


Nov 25 


53334.06 


138.1 


0.884-2.406 


WHT 4.2m 


2005 


Apr 19 


53479.86 


283.9 


0.477-0.944 


WHT 4.2m 


2005 


May 12 


53502.15 


306.1 


0.827-2.400 


Hiltner 2.4m 


2005 


Oct 14 


53657.08 


461.1 


0.889-2.453 


WHT 4.2m 


2005 


Oct 20 


53663.21 


467.2 


0.417-0.928 


WHT 4.2m 


2006 


Jan 15 


53750.03 


554.0 


0.890-2.343 


WHT 4.2m 


2006 


Dec 22 


54091.5 


895.5 


0.458-0.723 


Keck 10m 


2006 


Jan 21 


54121.5 


925.5 


0.444-0.890 


Keck 10m 



Note. — 

* We adopt an explosion date of 2004 Jul. 10.0 (MJD=53196.0). 

Column 5 details; 

DDO 1.88m: 1.88m telescope of the David Dunlap Observatory fVi nko et ani200g i 

Lick 3m: Shane 3m telescope at Lick Observatory (Leonard ct al. 20 061) 

LAT Scorpio: Scorpio on the 6m large azimuthal telescope tChugai et al.ll2005l : IChugaill2006l ) 

WHT 4.2m; 4.2m William Herschel Telescope of the Observatorio del Roque de los Muchachos. Slit width = 1". 
Hiltner 2.4m: 2.4m Hiltner Telescope of the MDM Observatory. 
Keck 10m; 10m telescope of the W. M. Keck Observatory. 



TABLE 5 

Optical and Near-Infrared Line Fluxes, F, and Luminosities, L. 





Ha 




Pa/3 


[O I] 6300 A 


[Fe II] 7155 A 


[Fe II] 12570 A 




Epoch (d) 


F 


L 


F 


L 


F 


L 


F 


L 


F 


L 


^rad 


89 


5300* (260) 


63(3) 


















586.9 


128 


2120* (100) 


25(1) 


















417.2 


138 






380(20) 


4.5(2) 














381.4 


283 


1535(80) 


18.1(9) 






250(10) 


3.0(2) 


82(14) 


1.0(2) 






101.7 


306 






137(7) 


1.62(8) 










35(10) 


0.4(1) 


81.5 


461 






39(4) 


0.46(5) 










20(2) 


0.24(2) 


16.4 


467 


230(10) 


2.75(15) 






134(7) 


1.58(8) 


41(7) 


0.49(8) 






15.35 


554 






9.7(8) 


0.11(1) 










12(3) 


0.14(3) 


6.0 


895 


10.7(7) 


0.13(1) 






7.6(9) 


0.09(1) 


3.0(2) 


0.035(3) 






0.18 


925 


9.5(7) 


0.12(1) 






6.9(3) 


0.082(4) 


2.52(8) 


0.031(1) 






0.14 



Note. — Line fluxes, F axe in units of 10~^^ erg cm^'^s"^ while dereddened line luminosities, L, are in units of 10'^* ergs~^. Figures in 
brackets give the error on the last one or more significant figures. The errors in the fluxes and luminosities are primarily due to uncertainties 
in the absolute fluxing and the levels of the underlying continua. Not included in the luminosity errors are systematic uncertainties in the 
distance and extinct ion (se e jjl.l) . 
* Spectra from lVfnko et al.l ((2009 ). 

In the final column is show n I/rad; the rad io active deposition pow er corresponding to the ejection of 0.0095 Mq of ^^Ni, scaled from the 
SN 1987A case specified bv lLi et al.l |[T9M ) fc lTimmes et al.l ifl^M ). 
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TABLE 6 

Optical and Near-Infrared Line Profile Parameters 



Epoch 




Velocity (kms ^) 




-Rbhm 


Model 


(d) 


'^BHM 


^pcak 


'^RHM ""HWHM 


(10^5 cm) 


(kms-l) 


Ha and Pa^ 


89* 


-2010 


-450 


1630 


1820 


1.5 




95t 


-2120 


-270 


2210 


2160 


1.7 




100** 


-2550 


-1620 


2000 


2270 


2.2 




128* 


-2670 


-1770 


-700 


980 


2.9 




138(Afi'ii;) 


-2190 


-1780 


-310 


940 


2.6 




283 


-1710 


-730 


990 


1350 


4.2 




306(Af/R) 


-1840 


-350 


910 


1370 


4.9 




461{NIR) 


-1170 


-180 


300 


730 


4.7 


-2400 


467 


-1240 


-160 


550 


890 


5.0 


-2400 


554(Af/_R) 


-800 


15 


500 


650 


3.8 


-2400 


895 


-310 


-140, 160 


-10" 


150" 


2.4 


-2400 


925 


-430 


- 80, 210 


_b 




3.4 


-2400 


[O I] 6300 A 


283 


-1760 


-140 






4.3 




467 


-1040 


-300 


400 


720 


4.5 


-2250 


895 


-440 


-210, 170 


-90" 


170" 


3.4 


-990 


925 


-440 


-170, 210 


20" 


230" 


3.5 


-990 


[Fe II] 7155 A & 12570 A 


283 


-1600 


-830 


370 


980 


3.9 




3D6(Af/ii;) 


-1670 


-510 


400 


1030 


4.4 




461{NIR) 


-1610 


-600 


160 


880 


6.4 


-3000 


467 


-1130 


-330 


510 


820 


4.6 


-2250 


554(Af/_R) 


-940 


-230 


260 


600 


4.5 


-2250 


895 


-400 


-210, 490 


-120" 


140" 


3.1 


-780 


925 


-410 


-170, 480 


-30" 


190" 


3.3 


-900 



Note. — All the velocities are with respect to the centcr-of-mass rest frame of SN 2004dj. 
Velocity subscripts: 
BHM: Blue wing, half-maximum 
RHM: Red wing, half-maximum. 
HWHM: Half-width, half-maximum. 

_Rbhm is the product of vbum and the elapsed time (epoch). 

*: IVinko et all i200^ . t : [Leonard eFaTI l l25ol) . ** IChugai eTaLl l l2005f) 

(NIR): Indicates NIR lines Pa/3 and [Fe II] 12570 A. 

In col. 3, 895 d and 925 d, the velocities of the stronger, blueshifted peak and the weaker, redshifted peak are given (see §2.4 and §3.4.4.1). 
": Based on blueshifted component only. 

Blueshifted component insufficiently resolved. 
Typical uncertainties in the period 89-554 days, in kms^^, are: ±50 (BHM), ±50 (peak), ±100 (RHM), ±60 (HWHM). At 895 d and 
925 d the uncertainties are smaller by a factor of ~ 4. 

The last column, "Model", gives the maximum velocity derived from profile (blue wing) matches (5 = 0)- see §3.4.4.1. 
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TABLE 7 

Parameters for Triple-Blackbody Matches to SN 2004dj Continda 



Epoch 


^hot 

(km s^-"-) 


^ hot 


^warm 

(kms^-"-) 


^warm 


warm 


^cold 

(kms^-*-) 


T 1 1 

^ cold 


-^hot 

CI 0^8 

prpr o 1 


warm 
PTP" R" 1 


r. rt 

-^cold -^total 
(1038 (1038 

ergs—-"-) ergs—-"-) 


-^rad 
('1038') 

orp" s ^ 


89 


450 


7000 


1750 


1.3 


1800 


7000 


200 


2141 


135 


0.33 


2276 


592 


106 


220 


7000 


1300 


1.2 


1800 


8500 


200 


695 


106 


0.69 


801 


507 


129 


180 


6000 


1250 


1.4 


1300 


7500 


200 


372 


39.5 


0.80 


412 


411 


251 


50 


6000 


580 


1.3 


1000 


5000 


200 


109 


11.3 


1.3 


120 


137 


281 


50 


5300 


470 


1.1 


900 


4500 


200 


82.9 


6.1 


1.4 


89 


103 


500 


3.8 


10000 


570 


2.5 


840 


1800 


200 


19.2 


21.5 


0.69 


41 


10.7 


652 


(1.0) 


10000" 


550 


3.1 


700 


1400 


200 


2.3 


16.4 


0.71 


19 


2.1 


859 


(0.16) 


10000" 


450 


3.3 


620 


1250 


200 


0.10 


11.7 


0.98 


12 


0.26 


996 


(0.05) 


10000" 


440 


3.8 


530 


800 


200 


0.013 


8.1 


0.5 


8.1 


0.077 


1207 


(0.01) 


10000" 


330 


3.4 


480 


1050 


200 


7.8 X 10-* 


4.5 


1.4 


4.5 


0.019 


1393 


(0.0025) 


10000" 


320 


3.8 


440 


900 


200 


6.5 X 10-5 


4.0 


1.3 


4.0 


0.009 


Note. 


— * The continuum luminosities tend to underestimate the total SN luminosity as 


they do not include line 


emission. 


This is 



particularly the case around 200 — 500 d when the relative contribution of nebular line emission is at a maximum. There is less of a problem 
before this era when the hot continuum dominates, or afterwards when the warm/cold continuum increasingly dominates. 
^- ^total excludes Lcold- This is because the cold component is due to an IS IR echo (see §3.4.3.1) which was predominantly powered by 
the peak luminosity of the SN prior to the earliest epoch of observation. 

**: Hot blackbody velocities and temperatures estimated by fixing the temperature at the 500 d value and extrapolating the earlier velocity 
evolution. 

In col. 13, I / rad th e radioactive depositi on corresponding to the ejection of 0.0095 Mq of ^^Ni, scaled from the SN 1987A case specified 
bv lLi et al.1 1(1993 1 fc lTimmes eTaLl lfl9M l. 

The 89 d and 996 d warm and cold parameters are based on photometry only. 



TABLE 8 

Hot Blackbody, Free-Free/Free-Bound and Warm Dust Disk (5 = 0) Parameters for Matches to SN 2004dj SEDs. 



Epoch 


'Jhotbb 


^hot 


''fb-ff 


Tfb-ft- 


»^c{fb-ff) 


dmax 


^d* 


T0.55 Mm 




T24 Mm 


A^dust 


(d) 


(km s" -"-) 


(K) 


(kms--"-) 


(K) 


(lO'' cm-3) 


(10^5 cm) 


(K) 








(IO-^Mq) 


89 


450 


7000 


3420 


(6000) 


100 














106 


275 


6000 


3420 


(6000) 


86 














129 


180 


6000 


3420 


(6000) 


56 














251 


50 


6000 


3080 


5000 


36 














281 


50 


5300 


2980 


5000 


24 














500 


9.0 


6000 


2470 


5000 


1.8 


5.2(5) 


650(15) 


9.6 


0.45(5) 


0.16 


0.22(2) 


652 


(1.0) 


(10000) 


2090 


5000 


0.6 


4.9(5) 


610(15) 


15.3 


0.73(10) 


0.26 


0.32(5) 


859 


(0.16) 


(10000) 


1590 


5000 


0.25 


4.6(5) 


570(15) 


17.4 


0.82(10) 


0.29 


0.33(5) 


996 


(0.05) 


(10000) 


1400 


2500 


0.14 


4(1) 


520(30) 


36.7 


1.7(5) 


0.62 


0.5(1) 


1207 


(0.01) 


(10000) 


1400 


2500 


0.05 


3.8(6) 


460(30) 


>80 


>4 


>1.5 


>1.0 


1393 


(0.0025) 


(10000) 


1400 


2000 


0.024 


3.8(6) 


430(20) 


>120 


>6 


>2 


>1.5 


Note. — 


"^dmax ' 


and Td are 


, respectively, the radius and temperature of the 


warm, dusty disk model. Uncertainties 


in the last one 



two figures are shown in brackets. Likewise for the estimated disk dust mass, M^^^^ (last col.). An example of uncertainties in the optical 
depths is shown in col. 10 for 10 /im. 

For the hot blackbody at 652 d and later, the temperature was fixed at 10,000 K and the velocity obtained by extrapolation, and so these 
parameters are shown in brackets. 

For the fb-ff modelling, the hydrogen temperature and electron density (cols. 5 & 6) for 251-1393 d were derived from the observed hydrogen 

velocities for SN 2004dj (col. 4) in conjunction with the late-time SN 1987A study of[Kozma & Fransson (1998) (see §3.4.3.2). For 89-129 d, 

the hydrogen temperature is shown in brackets as it is a rough estimate (see §3.4.3.2). 

Cols. 9-11 give the optical depths perpendicularly through the dust disk at the wavelengths indicated. 

The matches also included a CDS IR echo at epochs 89 — 500 d and an IS IR echo at all epochs (see §3.4.2 and §3.4.3.1). 
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TABLE 9 

Hot Blackbody, fb+ff and Warm (5 = 0) Dust-disk Luminosities compared with Radioactive Input. 

Fixed disk radius 



Epoch Li-jot -J^CDS ifl+fb ^IDDM ^total -^^rad 

(d) (10^8 (lO^s (lO^s (lO^s (lO^s (lO^S) 

ergs^-*^) ergs""*^) ergs"-'^) ergs"-"^) ergs~^) ergs"-*^ 



89 


2050 


128 


8.0 




2186 


592 


106 


586 


55.7 


10.0 




652 


507 


129 


372 


26.8 


7.7 




407 


411 


251 


109 


6.6 


3.3 




119 


137 


281 


82.9 


4.8 


4.2 




91.9 


103 


500 


14.0 


0.6 


0.20 


12.4 


27.2 


10.7 


652 


(2.3) 




0.03 


9.9 


12.2 


2.1 


859 


(0.10) 




0.005 


6.9 


7.0 


0.26 


996 


(0.032) 




0.0025 


3.9 


3.9 


0.077 


1207 


(7.7 X 10-*) 




0.0006 


2.3 


2.3 


0.019 


1393 


(6.5 X U)-'') 




0.0002 


1.8 


1.8 


0.009 



Note. — * Liddm is the luminosity of the warm, dusty disk modeL 
For the hot blackbody (col. 2) at 652 d and later, the temperature was fixed at 10,000 K and the velocity obtained by extrapolation, and 
so the luminosities derived from these parameters are shown in brackets. The post-652 d contribution of the hot blackbody to the total 
luminosity is negligible. 

In col. 7, Lrad is the radioactive deposition corresponding to the ejection of 0.0095 Mq of ^^Ni, scaled from the SN 1987A case described 
bv lLTet al. (1993) & Timmcs ct al. (1996). 

No IS IR echo luminosities are shown since these were predominantly powered by the peak luminosity of the SN prior to the earliest epoch 
of observation. In the wavelength region (A > 15 /im) where the cold component makes a significant contribution (> 20%) to the total 
flux, the IS IR echo model described in the text maintained a near-constant luminosity (15—150 fim) of ~ 2.0 X 10^® erg s~*between 89 d 
and 1393 d. 
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Fig. 1. — Sequence of images at 8 fim at 257 d (LHS), 621 d (middle) and 996 d (RHS). SN 2004dj is the point source at tlie center of 
eacli image. The supernova is clearly brighter on 621 d. This is due to the epoch being close to the peak of the thermal emission from the 
ejecta dust at this wavelength (see Fig. [J}. The fields are about 2 arcmin across (~ 1.8 kpc at the distance of the supernova). North is 
~ 30° clockwise from the upward vertical 




200 400 600 800 1000 1200 1400 

time after explosion (days) 



Fig. 2. — MIR light curves of SN 2004dj. They are uncorrected for reddening. For clarity, the plots have been shifted vertically by the 
amounts shown in brackets (mjy). Off to the right of the plot there are four additional 3.6 /im points spanning 1954 d to 2143 d (Table[TJ. 
They show little significant change during this period, with a mean flux of 0.32 it 0.02 mJy. This is consistent with there having been no 
decline since 1372 d. The MIR, flux of S96 alone has not been measured and so the light curves are uncorrected for 896. An estimate of 
the S96 contribution is given in Table [T] and is discussed in §3.1. 
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Fig. 3. — Evolution of the MIR spectral energy distribution of SN 2004dj. The SEDs are uncorrected for reddening. For clarity, the 
individual SEDs have been scaled by the amounts shown in brackets. As in Fig. [2] the fluxes from 896 have not been subtracted. The large 
peak at 4.5 ^tm present at 114 d and 257 d is due to the dominance of CO fundamental emission in this band during this period. 
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Fig. 4. — Optical and NIR light curves of SN 2004dj. For clarity, they have been displaced vertically by the values shown in brackets. 
The optical light curves are from Vinko ct al. (2006, 2009) (solid dots). In addition the V and R points of Zhang et al- (2006) (open circles) 
around the end of the plateau have been added to enhance the detail of this phase. For all data sets the epoch has been adjusted to our 
explosion epoch of MJD=53196.0. The NIR data are from the present work. Also shown is the 3.6 light curve from the present work. 
The vertical dashed line is at the epoch of the earliest 3.6 fim point and indicates the corresponding phase in the other light curves. S96 
fluxes at BVRI and at JHKs have been subtracte d from th e data. Also s hown for compa rison (labelled "Rad.") is the temporal evolution 
of the radioactive energy deposition as specified bv lLi et al.l I11993i ) fc iTimmes et all 119961) for SN 1987A with the addition of the early-time 
contribution of ^^Ni assuming complete absorption. 
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Fig. 5. — Evolution of the MIR spectra of SN 2004dj. The spectra have not been dereddened and the S96 flux has not been been 
subtracted. No correction has been apphcd for the hehocentric velocity of the SN. The 1207 d spectrum has a small exposure time and 
so is of relatively low S/N compared with those of 859 d and 1393 d (see Table [Sj. Also shown are the fluxes obtained through aperture 
photometry of the IRAC, PUI and MIPS images, interpolated to the epochs of the spectra. These photometric data are represented by 
round or square points respectively on alternate epochs. All plots except the latest have been displaced vertically for clarity, with the zero 
flux levels indicated by the horizontal dotted lines on the left. The stronger features are identifled, with the flducials redshifted to the SN 
rest frame. The wavelengths are labelled in microns. 
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Fig. 6. — Evolution of t he optical spectra of SN 2004dj. The 128 d spectrum is from lVinko et al.l H2006l 'l and the 895 d spectrum has 
already been presented in Vi nko et al.l f200S). The spectra have not been dereddened. In addition, the S96 flux has not been subtracted 
and this accounts for the relatively large continuum fluxes at epochs 895 d and 925 d. No correction has been applied for the heliocentric 
velocity of the SN. All plots except the latest have been displaced vertically for clarity, with the zero flux levels indicated by the horizontal 
dotted lines on the left. Note that the zero flux shown at 23 mjy corresponds to the 895 d spectrum. Also, the spectra have been flux-scaled 
by the amounts shown in brackets. Locations of spectral lines of interest are identifled, with the flducials redshifted to the SN rest frame. 
The wavelengths are labelled in microns. 
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Fig. 7. — Evolution of the NIR spectra of SN 2004dj. The spectra have not been dereddened and the S96 flux has not been subtracted. 
No correction has been applied for the heliocentric velocity of the SN. All plots except the latest have been displaced vertically for clarity, 
with the zero flux levels indicated by the horizontal dotted lines on the left. Also, the spectra have been fiux-scalcd by the amounts shown 
in brackets. The stronger features are identified, with the fiducials redshifted to the SN rest frame. The wavelengths are labelled in microns. 
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Fig. 8. — Evolution of spectral line luminosities of SN 2004dj. The luminosities h ave been correc ted for extinction. Also shown is the 
deposition power of 0.0095 Mq of 56Ni, scaled from the SN 1987A case described bv lLi et all 1 11993 ) &|Timmcs ct al. (1996). The errors 
in the fluxes and luminosities are primarily due to uncertainties in the absolute fluxing and the levels of the underlying continua. Not 
included in the luminosity errors are systematic uncertainties in the distance and extinction (see §1.1). 
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Fig. 9. — Evolution of spectral line velocities of SN 2004dj (shifted to the center-of-mass rest frame of the SN). The plots are as follows: 
triangles=red wings, half-width half-maximum, round dots=peak intensity, squares=blue wings, half-width half-maximum. On 895 d and 
925 d, the red wing half-width half-maximum refers to the blueshifted component only. Typical uncertainties in the period 89-554 days, 
in kms-l, are: ±100 (RHM), ±50 (peak), ±50 (BHM). At 895 d and 925 d the uncertainties are smaller by a factor of ~ 4. 
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Fig. 10.— Evolution of the spectral li ne profiles of SN 2004 di. The 89 d and 128 d s pectra are from lVinko et al.l I^OOB ). while the 95 d 
and 100 d spectra are from, respectively. [Leonard et aU l|200g) and lChugai et al] 1 12003 ). All other spectra are from the present work. The 
horizontal axes are in terms of equivalent velocity with respect to the SN center of mass, which has a velocity of +221 kms~^ relative to the 
Earth. The LH panel shows Ha and Pa/3 profiles. The Pa/3 profiles are labelled "IR". The middle panel shows [O I] 6300 A ([O I] 6364 A 
can also be seen lying at ~+3000 kms-^in the [O I] 6300 A rest-frame plots). The RH panel shows [Fe II] 7155 A and [Fe II] 12567 A 
profiles. The [Fe II] 12567 A profiles are labelled "IR". All the profiles have been flux scaled and shifted vertically for clarity. 
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Fig. 11. — Two-blackbody representation of the S96 SED. The sohd and open dots represent the pre-cxplosion and post-explosion fiuxes 
respectively ( Vinko et al. 2009). The latter data were acquired at 800 d. In the model, the fiuxes longward of ~1 /im are dominated by the 
cold component, with a temperature of 3500 K. The temperature of the hotter component, 50,000 K, is not intended to have a particular 
physical meaning, but simply serves as a means of representing and extrapolating the optical SED. 
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Fig. 12. — BLCs of IVinko et al.l I I200 6I) (open circles). IZhang et al.| 1 120061) (open triangles) and the hot+warm luminosities obtained 
in t he present wo rk via blackbody continuum matching (solid squares) (see §3.3 and Table [T] col. 11.). The phases of the Vinko ct Ed] 
and IZhang et all BLCs have been shifted to an explosion date of MJD=53196.0, the luminosities scaled to 3.13 Mpc, and the e xtinct ion 
to Ay = 0.31 mag. These BLCs are compared with the radioactive deposition power for SN 1987A as specified by ILi et al] 1 11993) & 
ITimmes et al.l ljl999 ). Two scaled cases are shown viz. 0.0095 Mq (solid line) and 0.016 Mq (dashed line) of S^Ni. We also show the total 
radioactive luminosity (red) for the 0.0095 Mq case. In addition (dotted line) are the actual UV-augmented bolometric light curves of 
SN 1987A llPun et al.lll995l) scaled to the case of an initial 56 Ni mass of 0.0095 M0. The divergence between the two SN 1987A datasets 
at late times is ascribed to differences in the IR flux measurements between the two observatories at, respectively, CTIO (lower curve) and 
ESQ (upper curve). 
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Fig. 13. — Comparison of blackbody continua with the optical-NIR-MIR photometric fiuxes and spectra of SN 2004dj at six epochs 
spanning 89 — 500 d. The solid dots show contemporary photometry. The blue, green and red lines show, respectively, the hot, warm and 
cold blackbodies. The cyan line shows the total model spectrum. 
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Fig. 14. — Comparison of blackbody continua with the MIR photometric fluxes and spectra of SN 2004dj at five epochs spanning 
652 — 1393 d. The solid dots show contemporary photometry. The green and red lines show, respectively, the warm and cold blackbodies. 
The cyan line shows the total model spectrum. A hot component was also included in the matches but had a negligible effect on the 
warm-cold matches to the MIR continua. Consequently, and in order to show the MIR behavior in more detail, the optical/NIR region is 
not shown. Estimates of the hot component contributions are given in Table [7] 
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Fig. 15. — CDS IR echo model compared with the 3.6 jum and 8.0 jum MIR excess Ught curves of SN 2004dj (LH panel), and the 89, 106 
and 129 d MIR excess SEDs (RH panel). The grain growth began at 50 d with a growth timescale also of 50 d. Amorphous carbon grains 
were adopted of radius 0.2 fxm.. The final dust mass is 0.33 X 10-5 M0. 
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Fig. 16. — Illustration of the line-profile models. The model profiles are derived from a homologously-expanding, emitting gas sphere 
attenuated by an embedded, opaque dust distribution. A minor additional contribution to the model is provided by a thin, expanding con- 
centric ring of emitting gas oriented with the normal to the ring plane at 6° to the line-of-sight. The gas-dust configurations are illustrated 
in the LH panels which show sections through the SN system in the plane defined by the thin ring axis and the line-of-sight. The gas sphere 
is shaded green and the attenuating dust shaded magenta. The thin ring, viewed edge-on in the figure, is represented by two large cyan dots. 
The region of line emission occluded by the dust is shown by the black shading. The expansion velocities are indicated by the axes. The 
final model profiles are obtained by combination of the attenuated sphere and unattenuated ring spectra. The resulting spectrum is then 
smoothed to the spectral resolution. The RH panels show the resulting model line profiles (solid red), plotted in velocity space, compared 
with the observed spectrum (blue). Also shown (dotted red) is the line profile that would result in the absence of dust attenuation. In this 
illustration we compare the model profiles with the observed 895 d [O I] 6300 A line. In this case, the gas sphere emissivity declines as r~^-^. 

The upper panels illustrate the case where the dust is distributed as an opaque, concentric sphere. Two sizes of sphere are considered. 
Profiles a and b (upper RH panel) correspond, respectively, to the large and small dust spheres (upper LH panel). It can be seen that 
neither reproduce the observed line profile. The lower panels illustrate the case where the dust is composed of amorphous carbon and is 
distributed as a uniform density, nearly face-on concentric disk with t = 18 perpendicularly through the disk plane. This opti cal depth is 
the value required for the model to simultaneously match the line profile and the contemporaneous MIR continuum (see Fig. I19|l . While 
the disk is shown as being coplanar with the ring, for ease of computation the disk is taken to be exactly face-on. Shifting the disk tilt to 
match that of the ring would have had only a small effect on the model profile. No correction is made for minor edge effects at the disk's 
outer limit (radius). Also shown (lower RH panel) are the intrinsic line profiles from the attenuated gas sphere (green solid line) and the 
thin ring (cyan solid line). The thin ring emission is not affected by the dust. For clarity, both of these intrinsic profile plots have been 
scaled downward in fiux by the same amount relative to the total profile. It can be seen that the final dust-disk model profile provides an 
excellent match to the data. 
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Fig. 17. — Individual dust-disk model matches (red) to observed line profiles (blue) of SN 2004dj for the case where 5 = (uniform dust 
density). The velocities of the observed spectra have been shifted by -221 kms~^ to match the center-of-mass rest frame of the SN. The 
latest two epochs also include a small contribution from a ring of emission. Also, the velocity scale of these two epochs has been expanded 
to provide more detail. The matches incorporate dust disk masses obtained by interpolating to the profile epochs the values obtained from 
the MIR continuum modelling. Between 554 d and 925 d the dust mass increased from 0.25 X 10""* Mq to 0.4 X lO""* Mq. The radius of the 
dust disk decreased slightly from 5.3 X 10^^ cm to 4.6 X 10^^ cm. The disk half-thickness expanded from 0.45 X 10^^ cm to 0.90 X 10^^ cm. 
For each profile the optical depth had a value in the range 3 < r < 19 (i.e. the disk was always of high optical depth). 
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Fig. 18. — 89-281 d: Model continua comprising hot blackbody (solid blue), free-bound (dashed blue), free-free (dashed red), CDS IR 
echo (green), IS IR echo (solid red) and total flux (solid cyan). The observed spectra and contemporary photometry (dots) of SN 2004dj 
arc plotted in black. 

500 d: As for earlier epochs but with the addition of IR emission from a warm disk of dust (magenta). The dotted cyan line represents the 
model continuum when the IDDM is replaced with a blackbody adjusted to match the shortwave region of the MIR continuum. At this 
epoch the ff and fb fluxes are too weak to appear on the plot. 
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Fig. 19. — Model continua at epochs 652 — 1393 d, comprising IR flux from a warm disk of dust (magenta), IS IR echo (red) and total flux 
(solid cyan). The dotted cyan line represents the model continuum when the IDDM is replaced with a blackbody adjusted to match the 
shortwave region of the MIR continuum. By 1393 d the blackbody has fully merged with the IDDM; i.e., the model is effectively opaque 
over the whole observed wavelength range. The observed spectra and contemporary photometry (dots) of SN 2004dj are plotted in black. 
A hot blackbody plus free-bound and free-free continua were also included in the matches but these had a negligible effect on the MIR 
continua. Consequently, and in order to show the MIR behavior in more detail, the optical/NIR region is not shown. Estimates of the hot 
and fb-ff component contributions are given in Tables |8] and |9] 



